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Cryogenic Propellant Scavenging 
Final Report 

August 1982 - March 1985 

B. Louie, N. J. Kemp, and D. E. Daney 

Chemical Engineering Science Division 
National Bureau of Standards 
Boulder, Colorado 80303 

This report is a detailed description of a computer model that has been 
developed for assessing the feasibility of low g cryogen propellant scavenging 
from the Space Shuttle External Tank (ET). Either pump-assisted or pressure- 
induced propellant transfer may be selected. The receiver tank is chilled by 
emitting a low flowrate of single-phase cryogen through small nozzles. When 
two phases are present the flowrate is increased to represent transfer through 
the main piping. The program will accept a wide range of input variables, 
including the fuel to be transferred (LOX or LH 2 ), heat leaks, tank 
temperatures, and piping and equipment specifications. 

The model has been parametrically analyzed to determine initial design 
specification for the system. Pressure-induced transfer of LH 2 can be 
accomplished in approximately 7 minutes with a 3-inch (0.076 m) line size. 
Pump-assisted scavenging of LOX can be completed in less than 4 minutes by 
using a 2 HP (1491.4 W) pump and a 3-inch (0.076 m) line. To maximize the 
quantity of LH 2 recovered, the receiver tank should be prechilled to -290 °F 
(94.1 K). It was determined that the LOX receiver tank does not require 
prechilling and can have a temperature as warm as 300 °F (421.9 K) without 
significant venting of fluid. 


Key words: computer model; cryogenic; fluid transfer; low-g; scavenging; 

thermodynamics. 
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1. INTRODUCTION 


The continuing success of the frequent Space Shuttle flights has encour- 
aged interest in suborbital scavenging of propellants. Currently, the surplus 
liquid hydrogen (LH^ and oxygen (LOX) are jettisoned with the external tank 
(ET) or vented to space after the main-engine cutoff (MECO). Nearly 10,000 
pounds* (4535.9 kg) of fuel remain after MECO, which can be compared to the 
full payload of 65,000 pounds (29483.4 kg). Significant savings can be 
realized immediately with the proposed scavenging. Accomplishment of this 
task would support the longer range concepts of tethered vessel storage or 
orbiting space station operations. 

Brux and Stefan [1] have studied the physical phenomenon of cryogenic 

propellant transfer and discuss the logistics of the task, including receiver 

tank configuration and location, various heat fluxes, and the duration of a 

mated coast period after MECO required to meet settling thrust requirements 

for the ET. They find the concept feasible, provided a positive acceleration 
-4 

of 1 x 10 g's or greater is continued for a period of five to twenty 

minutes. A period of twenty minutes is required for the worst-case transfer 
of both fluids. The transfer of liquid oxygen requires a 2 to 7. 5- inch (5 to 
19 cm) line size and a boost pump. Liquid hydrogen transfer can be 
accomplished using a 4-inch (10 cm) line size with no pump. Simulation of the 
transfer incorporates a receiver tank cooldown period, during which a low 
flowrate of cryogen is circulated. 

*This report does not conform with NBS policy on SI units. At the request of 
the sponsoring agency, NASA-JSC, English units are used for the input and 
output of data. 
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The present investigation was undertaken to determine the possibility of 
propellant scavenging and to evaluate the conditions required for successful 
ET residual transfer. A computer model written in FORTRAN IV (FORTRAN 66) is 
based on a thorough thermodynamic analysis of the shuttle receiver vessel, 
including cooldown of the vessel interior, and heat fluxes to the ET and 
shuttle tank. The model is analyzed parametrically in order to provide 
information to be used in the design of such a system. 

This report summarizes our development of a computer model for cryogenic 
propellant transfer at the National Bureau of Standards. Fundamental 
equations, system description, and assumptions are presented in section 2 
entitled PROGRAM DEVELOPMENT. Program logic and descriptions of parameters, 
subprograms, and numerical methods are discussed in section 3 entitled PROGRAM 
DESCRIPTION. The effect on the overall time of transfer and other operational 
characteristics is analyzed using parametric variations. Factors that are 
considered include heat leaks, vessel prechi 11 ing, and cooldown requirements. 
These topics are discussed in section 4 entitled RESULTS AND DISCUSSION. 
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2. PROGRAM DEVELOPMENT 


The physical situation described by the computer model is the transfer of 
cryogenic propellant between two tanks in a low-gravity environment. A simpli- 
fied schematic of the propellant scavenging system is shown in Figure 1. The 
points labeled 1 and 2 respectively denote the main fill line into and vent 
line from the receiver tank. 



Figure 1. Flow schematic for cryogenic propellant scavenging 
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The process begins with a cooldown of the receiver tank by spraying the tank 
walls with cryogen. This is accomplished by emitting a low flowrate of 
cryogen through tangential spray nozzles while the main fill line is closed. 
Following the cooldown period the cryogen is transferred through the main fill 
line by utilizing ullage pressure in the supply tank or by pumping. The 
transfer is complete when the receiver tank is full, when the supply tank is 
empty, or when the back pressure in the receiver tank becomes too large. 

The model is developed from the differential mass and energy balances 
applied to the receiver and external (ET) tanks. There are two separate 
thermodynamic conditions which must be considered for the receiver tank. 

During the cooldown period, there is only vapor in the tank. Following this 
single-phase process, a two-phase transfer into the tank by the cryogen must 
be described. Allowances for pressure relief are made to cover those 
instances where the tank pressure becomes too great and venting of the vapor 
to space results. The receiver tank derivation is based on the following 
assumptions: 

(1) Heat leaks into receiver and supply tanks (dQ and dQ $ ) are inde- 
pendent of tank temperature and pressure. 

(2) The enthalpy of the vented fluid is equal to the receiver tank 

vapor enthalpy (h v ). 

(3) Receiver tank vent pressure is independent of tank temperature. 

(4) There is sufficient settling so that only liquid is drained from the 
supply tank, and only vapor is vented from the receiver tank. 

A complete derivation of the single-phase receiver tank equations is 
given in Appendix A. Although the equations are derived for a differential 
mass, dm, they are equally valid when dm is considered a flowrate. Figure 2 
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shows the thermodynamic properties for the receiver tank associated with the 
fluid in the inlet, vent and tank. The equations are summarized as follows 
(definitions for the equation variables are found in the table of 
nomenclature): 

Mass balance - 



Figure 2. Flow schematic for the receiver tank 
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A complete derivation of the two-phase receiver tank equations is given 
in Appendix B. These equations are summarized as follows: 

Mass balance - 

dm^ + dm 2 = dm y + dm £ (7) 

Energy balance - 

dU = dQ + h-^dm^ + h^dm^ (8) 
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Constant volume - 


V = V 2 + \ "v 
0 = v^dm^ + m^dv^ + v y dm y + m v d\>. 


The density, p 0 or p , is substituted for the specific volume. Solving for 

X V 

the change in the mass of the liquid yields the following equation: 


CL 

(dm + dm ) - p 

[m ^ ♦»,, 

^ ldP 

?v 

1 2 £ 

2 3 p , 

8P V. 


If the pressure is less than the vent pressure, the change in the pressure is 
given by: 

f Po 1 / P 0 \ 
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dm 2 = 0. 
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If the pressure is the same as or exceeds the vent pressure, the tank 
pressure is constant. The differential pressure and mass flow through 
the vent are as follows: 


dP = 0 


u 


dm 2 = 


v 


Z xvSL 


- u v + V 1 - # 


dm^ + 


dQ (1 


\>Z J 


(h v - V T 


z 


( 12 ) 


The ET model is based on the assumption of non-equilibrium (thermal 
stratification) between the liquid and vapor phases remaining in the supply 
tank. This assumption is valid because the prechilled supply tanks contain 
subcooled liquids and the nominal amounts of fuel remaining after MECO are 
extremely small. Thus, no significant mass transfer can take place between 
the phases unless the temperature increases enough to vaporize the liquid. 
When this occurs, the transfer of fuel stops since severe pump cavitation 
conditions are present. 

The following set of equations describes the supply tank at any time: 


Energy balance - 
vapor phase: 

liquid phase: 
Pressure equation 


dq v 

d \ = if + v p 


dh * - i-j* + v p 


dl V , V q v , 


dQ, 


dP = 


Vv 6 v 


9 A 


v 


£ + Vv 2 


e A "aWv 


(13a) 

(13b) 


(14) 
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Temperature equation - 




» A dP 


(15) 


Derivations of the above equations can be found in Appendix E. 

Once the supply tank liquid enthalpy is determined, the enthalpy of the 
fluid entering the receiver tank can be calculated by: 


h - h + ^i + w. 

"l s dn^ dm 1 


(16) 


and 

dh = dQ 
s x s 
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3. PROGRAM DESCRIPTION 


3. 1 Overview of the Computer Model 


Figure 3 shows the components of the computer simulation program and the 
communication paths as indicated by the direction of the arrows. 



Figure 3. Components of the computer model 

The program is written in FORTRAN IV and uses program subsections called 
subroutines and functions. The four major program sections are the main 
program SCAVAG and the subroutines COOL, SUPPLY and PFAIL. SCAVAG performs 
three distinct operations in modeling the scavenging of cryogenic propellant: 
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1) it serves to initialize the system parameters and direct the program flow; 

2) it transfers the calculations of the single-phase receiver tank and the 
external tank to subroutines; and 3) it performs the calculations of the 
two-phase receiver tank. COOL and SUPPLY perform the thermodynamic 
state calculations of the single-phase receiver tank and external tank, 
respectively. PFAIL sets up plotting arrays and prints messages signifying 
which completion procedure is chosen in the program. 

Provision is made at the start of the main program to select the propel- 
lant for the simulation. The namelist option is used to select the 
propellant. The program default value is NGAS=1, which signifies that LH 2 has 
been chosen. By using the namelist FUEL and setting NGAS=2, LOX can be 
selected. 

The default values for the parameters for the entire system are initial- 
ized according to the propellant chosen. Many of these values are the same as 
those used by Brux and Stefan [1]. The parameters and their default values 
are given in Table 1. A namelist option called NAME is used to change any of 
the initialized values. 

COMMON statements are used to pass the same value for a parameter between 
different portions of the program. For instance, the calculations in the main 
program and subroutines are dependent on the elapsed time or the quantity of 
propellant. Table 2 gives a general description of the parameters associated 
with each of the COMMON statements. 

Propellant selection is important in determining the correct fluid 
thermodynamic properties. The properties data are part of the Fluids Pack 
computer programs [2] which use routines based on real fluid properties for 
oxygen and parahydrogen to calculate various thermodynamic properties [3, 4, 

5]. The thermodynamic properties, such as saturation pressure, density, 
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Table 1. System Parameter Descriptions and Initialization Values* 


Parameter Value 


Location 

Name 

Description 

LOX 

lh 2 


Supply 

PET 

Pressure 

20 (0.138) 

32 

(0.221) 

Tank 

TET 

Temperature 

-315 (80.2) 

-425 

(19.1) 


QET 

Heat leak 

0 (0) 

90000 

(26376.3) 


VET 

Volume 

19786 (560260.4) 

53518 

(1515415.7) 


TOTAL 

Mass of propellant 

6270 (2844.0) 

3098 

(1405.2) 



to be transferred 




Receiver 

PRESS 

Pressure 

1 (0.007) 

18 

(0.124) 

Tank 

TEMP 

Temperature 

60 (288.6) 

-290 

(94.1) 


Q 

Heat leak 

0 (0) 

0 

(0) 


VOL 

Volume 

300 (8494.8) 

780 

(22086.5) 


MWALL 

Mass of tank wall 

350 (158.8) 

975 

(442.4) 


PVENT 

Vent pressure 

30 (0.207) 

30 

(0.207) 

Transfer 

DIAM 

Pipe diameter 

4 (10.16) 

5 

(12.70) 

Piping 

LENGTH 

Length of straight pipe 

100 (30.48) 

70 

(21.34) 


NELBOW 

Number of elbow fittings 

20 

10 



NGATE 

Number of gate valves 

0 

0 



NGLOBE 

Number of globe valves 

2 

2 



NANGLE 

Number of angle valves 

1 

0 



NBUTT 

Number of butterly valves 

0 

0 



METER 

Number of flow meters 

1 

1 



QPIPE 

Heat leak 

104400 (30596.5) 

50400 

(14770.7) 


POWER 

Pump power 

4 (2982.8) 

0 

(0) 

Cooldown 

CDLENG 

Length of straight pipe 

20 (6.10) 

20 

(6.10) 

Piping 

NNOZ 

Number of nozzles 

6 

6 



NOZDIA 

Nozzle diameter 

0.25 (6.4xl0- 3 ) 

0.500 

(1.3xl0- 2 ) 


HEADIA 

Header diameter 

1.000 (2.5xl0- 2 ) 

2.000 

(5.1xl0- 2 ) 


CONST 

1. /orifice constant 

1.639 

1.639 



*Units used to initialize parameters: pressure - psia (MPa), temperature - 
°F (K) , heat leak - Btu/hr (W), volume - cubic feet (liters), mass - pounds 
(kg), power - horsepower (W), length - feet (m), diameter - inches (cm). 
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and thermal conductivity, are found by using subroutines and functions built 
into Fluids Pack. When LOX or LH 2 is chosen, the subroutines DATA 02 or DATA 
PH2 are called. 

The distinction between propellants is considered at other points in the 
computer code. The plotting and numerical integration routines are dependent 
on the propellant chosen. PFAIL writes the propellant name into a data file 
which is used in producing the output plots. The time step in each numerical 
integration is based on the values of the flowrate into the receiver tank and 
which propellant is used. Also, since the receiver tank construction 
materials are different for LOX and LH 2 , the heat capacity calculation for the 
tank wall in FUNCTION CWALL requires this distinction. 


Table 2. Description of Parameters in the COMMON Statements 


Name 

Description of parameters 

PARAM 

initialized, constant values 

SUB 

time dependent 

FLOW 

time dependent and transfer piping 

NOZZLE 

cool-down piping 

PUMP 

pump values 

NGAS 

choice of propellant 

ETANK 

supply tank 

PROB 

programmed exit designation 


After the parameters have been initialized, they are converted to SI units. 
This is done to ensure compatibility with the thermodynamic data required from 
the Fluids Pack routines. All subsequent calculations are performed in SI 
units, with the exception of those in the plotting routine. 
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Before the transfer of fluid begins, the equivalent length of pipe between 


the supply and receiver tanks is calculated. Data for equivalent lengths of 
various pipe fittings are taken from Bennett and Meyers [6]. The equivalent 
length of each fitting is linearly proportional to the pipe diameter as given by 



D + C 


II 


(17) 


The constants are shown in Table 3. The total equivalent length is calculated 
by summing all the individual equivalent lengths. The equivalent length is 
then used by FUNCTION FLO and FUNCTION PUMPFLO to calculate a flowrate between 
supply and receiver tanks. 


Table 3. Calculated Constants for Equivalent Length 



c i 

c n 

90° elbow 

2.53 

0.414 

Gate valve 

0.571 

-0.0125 

Globe valve 

38.1 

0.107 

Angle valve 

14.2 

-1.18 

Butterfly valve 

1.79 

-0.147 


After parameter initialization and unit conversion has been accomplished, 
the thermodynamic condition of the supply tank is established. This procedure 
involves determining the quantity of ullage vapor in the tank. The initial 
flowrate between the tanks is calculated using the pressure difference between 
the tanks and the pressure drop through the piping with FUNCTION FLO. If the 
flow is aided by a pump, the increased flowrate as calculated in FUNCTION 
PUMPFLO is used. 
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3 . 2 Description of Program Sections 


SCAVAG 

The main program SCAVAG directs the sequence of calculations in the simu- 
lation. Following parameter initialization and unit conversion, the cool-down 
of the receiver tank takes place by calling SUBROUTINE COOL. When the receiver 
tank is cooled, there are two phases present and a return is made to SCAVAG. 

It is assumed that the transfer piping is simultaneously chilled. 

The remaining propellant is transferred in two phases. First, the 
enthalpy of the fluid entering the receiver tank is computed by using an 
energy balance on the transfer piping. The density and enthalpy of the liquid 
and vapor are found with Fluids Pack functions. These values are used in the 
first law of thermodynamics to find the internal energy of both phases in the 
receiver tank. Next, the integration of the differential pressure and mass 
equations, Eq. (10)-(12), is carried out over an incremental time step 
(forward Euler method) and incorporates the thermodynamic values above. A 
check is made to see if the pressure in the receiver tank exceeds the vent 
pressure; the mass of propellant in the tank is adjusted accordingly. 

Following this, the supply tank properties are updated by calling SUBROUTINE 
SUPPLY. Finally, the values for the pressure, flowrate, and propellant vapor 
and liquid quantities are stored in arrays for plotting. 

This calculation procedure is repeated within a loop until the receiver 
tank is full, until all the propellant has been transferred, or until another 
condition has been encountered which causes the program to be terminated. For 
instance, should the pressure in the tank exceed that of the entering stream, 
the transfer stops. If the simulation proceeds normally, SCAVAG calls 
SUBROUTINE PFAIL to prepare the data for graphic display. 
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SUBROUTINE COOL 


The single-phase cool-down of the receiver tank is performed by SUBROUTINE 
COOL until the receiver is at the correct pressure and temperature for a two- 
phase propellant. Contained in the call and subroutine statements are three 
parameters that are needed to transfer into and return from COOL. SCAVAG 
provides values for the mass of vapor in the tank and the mass of vented vapor. 
COOL returns to SCAVAG the differential pressure and the current values of the 
mass of vapor in the tank and the quantity vented. A calculation procedure 
similar to that in SCAVAG is followed to solve Eq.(l)-(6). 

Several thermodynamic quantities are needed to solve the single-phase 
differential equations. COOL calls SUBROUTINE PPHI which calculates the 
Griineisen parameter for the single-phase (vapor) propellant in the receiver 
tank. A heat balance analysis on the transfer pipe affords the enthalpy for the 
entering stream. The Fluids Pack functions are used to find the density, 
enthalpy, and specific heat of the fluid in the receiver tank. When all of 
these thermodynamic values are known, the equations are solved incrementally, 
with the appropriate checks for venting and back pressure. If all the propellant 
has been transferred before any liquid is present in the receiver tank, a call 
is made to PFAIL to prepare the plotting data and terminate the program. Again, 
the supply tank fluid properties are updated by calling SUPPLY. 

To determine when two phases are present in the receiver tank, saturation 
and actual gas densities are compared. If the gas saturation density at the 
receiver tank pressure is greater than the actual tank density, the receiver 
tank is two-phase. SUBROUTINE COOL returns to the main program following this 
determination. 
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SUBROUTINE SUPPLY 


When updated supply tank conditions are required in either the main pro- 
gram or SUBROUTINE COOL, SUBROUTINE SUPPLY is accessed to perform the necessary 
thermodynamic calculations. The fraction of the supply tank occupied by either 
phase is determined first so that the heat flux to the tank can be propor- 
tioned accordingly. The vapor phase is assumed to be at the tank pressure and 
slightly above the saturation temperature, thereby ensuring that condensation 
will not occur. Next, the Fluids Pack routines are called upon to provide 
thermodynamic values for the vapor and liquid phases. The Griineisen parameter 
and heat of expulsion for both phases are calculated by using SUBROUTINE PPHI 
and FUNCTION TH, respectively. Finally, the incremental changes in the tank 
pressure and liquid temperature and enthalpy are found by using Equations 
(14), (15) and (13b) and the appropriate time step. 

SUBROUTINE PFAIL 


SUBROUTINE PFAIL, with ENTRY FULL and ENTRY DONE, terminates the program 
and sets up the data and parameters to be plotted. The different entry points 
are accessed by unique call statements. PFAIL is called when: 1) the flow 
between tanks has stopped due to back pressure from the receiver tank; 2) when 
receiver and supply tank pressures are equal and a pump is not used; or 3) 
when there is excessive pressure in the receiver tank and a pump is used. 

PFAIL writes a message indicating why the simulation stopped and what the 
final conditions are. The other entries into PFAIL address other conditions 
and bypass inappropriate operations. ENTRY FULL is called when the main pro- 
gram or SUBROUTINE COOL determines that the receiver tank is full of propel- 
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lant. When all the fluid has been emptied out of the supply tank, ENTRY DONE 
is called. ENTRY PDATA is called when a programmed exit has been executed. 

Each entry into the subroutine follows a similar computational procedure. 
Initially, the units for the parameters to be printed in the termination mes- 
sage are converted to English from SI. Following this conversion, the message 
is printed on the output file stating the cause for program termination, the 
quantity of propellant transferred, and the temperature and pressure for each 
tank. The quantity of fluid vented from the tank is also given. 

Preparation of the data and parameters for plotting involves several 
steps. The upper and lower limits for the plots are determined first. The 
descriptions for the plotting parameters are shown in Table 4. Since each of 
the plots depicts a parameter as a function of time, the lower value is set to 
zero and the upper value is either 10 or a multiple of 50 minutes. Should the 
total time required for transfer be more than 10 minutes, the greater value 
will be used. 

The lower limit for the receiver tank temperature is based on the initial 
temperature in the supply tank, and the upper limit depends on the initial 
receiver tank temperature. Both values are multiples of 100 and in units of 
Fahrenheit. 

The limits on the receiver tank pressure are zero and 50 psia (0.345 
MPa), unless a value in the pressure array exceeds the upper limit. When this 
is the case, the value is incremented by 50 psia (0.345 MPa) iteratively until 
the maximum pressure data point is less than the newly set limit. The upper 
limit for the transfer flowrate is 2000 lb/min (907.2 kg/min) for LH 2 or 
10,000 lb/min (4535.9 kg/min) for L0X. 
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Table 4. Variables Used to Plot System Parameters 


Simulation variable 

Units 

Array 

Location of low, 
high bounds 

TIME 

Mi nutes 

X 

0PL0T(1) , 0PL0T(2) 

TEMP 

°F 

Y1 

0PL0T(3) , 0PL0T(4) 

PRESS 

psia 

Y2 

QPL0T(5) , 0PL0T(6) 

FL0W1 

Pounds/mi n 

Y3 

0PL0T(7) , 0PL0T (8) 

MVAP 

Pounds 

Y4 

0PL0T(9) , 0PL0T(10) 

MLIQ 

Pounds 

Y5 

0PL0T(11) , 0PL0T(12) 

TOTAL 

Pounds 

Y6 

0PL0T(13) , 0PL0T(14) 

ETQ 

Dimensionless 

Y7 

0PL0T(15) , 0PL0T(16) 

TANK 

Dimensionless 

Y8 

0PL0T(17) , 0PL0T(18) 


The quantity of propellant in the receiver tank is graphically displayed 
by plotting the vapor, liquid, and combined amounts over time. All these 
limits are dependent on the initial quantity of propellant (TOTAL) to be trans- 
ferred. The upper limit for the vapor quantity is initially set to 500 lbs 
(226.8 kg), but it is increased iteratively by 500 lbs (226.8 kg) if one- tenth 
of TOTAL is larger than the limit. The upper limits for the liquid and 
combined quantities are 10 times the value determined for the vapor quantity. 
The quality of the propellant in the supply tank and entering the receiver 
tank is stored in the array as percent vapor. The upper limit for the quality 
is 100% vapor and the lower limit is 100 percent liquid. 


3.3 Descriptions of Other Subroutines 


SUBROUTINE PPHI 


The thermodynamic quantity called the Gruneisen parameter, 4> , appears 
frequently in single-phase flow [7]. It is found in the equations describing 
the single-phase cooldown of the receiver tank (Eq. (3)-(4)) and the subcooled 
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liquid propellant in the supply tank (Eq. (14)) in the present study. Both 
subroutines COOL and SUPPLY call PPHI during their computational procedures 
for values of <j>. 

The dimensionless Griineisen parameter is given by 


or 




(18) 


(19) 


When the second form of <() is used, the Fluids Pack subroutine DPDT can be used 
to find the partial derivative and the constant volume specific heat of the 
fluid. The density and temperature are passed from either COOL or SUPPLY to 
PPHI which returns a value for <|>. 


SUBROUTINE CHECK 


Calls to this subroutine are routinely made to determine the quality of 
propellant. Both SCAVAG and COOL call CHECK to find the quality of the fluid 
entering the receiver tank. SCAVAG, COOL, and SUPPLY utilize CHECK to find 
the amount of vapor in the supply tank fluid. Because pump cavitation can be 
severe with a small quantity of vapor, the run will be terminated with this 
determination. A message is printed which states the percentage of vapor in 
the fluid. 

For inputs of pressure, temperature and enthalpy, SUBROUTINE CHECK deter- 
mines if the fluid is single-phase or two-phase and the mole fraction in each 
phase. First, the liquid saturation enthalpy is obtained from Fluids Pack 
using the input pressure. Next, if the input enthalpy is less than liquid 
saturation enthalpy, the stream is subcooled, and CHECK returns an indication 
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that no vapor is present (i.e. x = 0). If the input enthalpy is greater than 
the liquid saturation enthalpy, the proportion of fuel in the liquid and vapor 
phases is determined by: 


h = h J,sat (1 ' x) + 




or rearranging, 


x ^£,sat h)/(h^,sat ^v,sat^ 


( 20 ) 


Since the flow equations were developed for an incompressible fluid, 
there may be significant error if the fraction of vapor is much greater than 
zero. 


3.4 Descriptions of Functions 

Several thermodynamic quantities in the form of partial derivatives are 
required in the two-phase calculations in SCAVAG. While Fluids Pack directly 
calculates some partial derivatives such as DPDT (which is used in determining 
<{>), other partial derivatives must be calculated using an iterative approxima- 
tion. The approximation is written so that thermodynamic properties from 
Fluids Pack may be utilized. 

FUNCTIONS DVDPV and DVDPL 


FUNCTION DVDPV and FUNCTION DVDPL calculate the partial derivative of 
specific volume with respect to pressure for saturated vapor and liquid re- 
spectively. The pressure is the required parameter in the function statement. 
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These partial derivatives are required by Eq. (10) for the two-phase integra- 
tion. The partial derivative is approximated by the equation: 


/3v\ ”sat< P + AP > ' ' , sat< P - AP > 

\^ P /sat “ 2AP 

The values for the specific volume at saturation conditions are the inverse of 
the density computed by Fluids Pack at the same temperature and pressure. The 
pressure increment AP is initially set at 0.05 atm (0.005 MPa) and decreased 
until the partial derivative changes insignificantly. 

FUNCTIONS DUDPV and DUDPL 

FUNCTION DUDPV and FUNCTION DUDPL, given the pressure as an input, deter- 
mine the partial derivative of internal energy with respect to pressure for 
saturated vapor and liquid, respectively, as required by Equation (lib). The 
partial derivative is approximated in a similar manner to DVDPV and DVDPL by: 


3u\ 

m 


u sat (P + AP) ' u sat (P ' AP) 


sat 


2AP 


( 22 ) 


The internal energy, u^, is calculated using enthalpy values from Fluids 
Pack at the saturation temperature and pressure and the (First Law) energy 
relationship as given by: 

u = h - (Pv) 


or (23) 

u = h - (P/p) 

Again, AP is initially 0.05 atm (0.005 MPa) and decreased until the desired 
accuracy is obtained. 
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FUNCTION DTDP 


This function determines the partial derivative of temperature with 
respect to pressure for the saturated fluid (liquid and vapor). The partial 
derivative is approximated by the expression 


9T\ _ T saA P ^ AP > - T sat< P ~ AP > 

9P/ sat 2AP 


(24) 


Pressure is required as an input for the function and is subsequently used in 
the Fluids Pack function FINDTV to obtain the saturated temperature. The 
temperature is found at a pressure that is a small increment above and below 
the input. The increment, AP, is decreased from an initial value of 0.05 atm 
(0.005 MPa) until the convergence criterium is met. 


FUNCTION ALPHAL 


The temperature equation in SUPPLY contains the term for the liquid phase 
bulk expansivity, a , which is defined as (7): 


a = - 
a si pvax 


(25) 


FUNCTION ALPHAL calculates this value by approximating the partial derivative 
in an iterative manner. This approximation is given by: 


(Bp\= p(T+AT) - p(T-AT) (26) 

' 3T/ p 2AT 

A temperature increment AT is used to find values of the density at tempera- 
tures slightly above or below the liquid temperature. The temperature incre- 
ment AT is initially set at 0.05 K and decreased until the partial derivative 
changes insignificantly from one iteration to the next. 
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FUNCTIONS FLO, FRICT, PUMPFLO, and PEFF 


The functions used to calculate the transfer flowrate are dependent on 
the pressure in the receiver tank and are themselves interdependent. These 
functions are used to calculate the flowrate of propellant between the supply 
and receiver tanks. FLO and PUMPFLO are accessed from SCAVAG and COOL. FRICT 
is used in FLO and PUMPFLO, and PEFF is used in PUMPFLO only. A description 
of each follows. 

FUNCTION FLO 

FUNCTION FLO calculates the flowrate between tanks for a given pressure 
drop. This function is used to initialize the flowrate, and the value it 
provides is an input for FUNCTION PUMPFLO. SCAVAG passes the receiver tank 
pressure through the function statement. 

The flowrate between the tanks is related to the pressure drop in the 
system. The overall pressure drop can be stated as the sum of the individual 
pressure drops as follows: 


AP 


system 


= AP m • (dm 2 ) 
main 


+ AP 


meter 


+ AP 


cooldown 


(dm 2 ) 


(27) 


where (dm 2 ) is the functional dependence, not the multiplier. The overall 
pressure drop is the difference between the pressures in ET and receiver tank. 
The pressure drop in the meter is five percent of the overall pressure drop. 
The remaining terms in Eq. (27) are functions of the flowrate. The pressure 
drop through the main transfer line and the cooldown line are given by: 


AP. 


main 


32Cf . L . (dm) 2 
mai n mai n 

pn 2 D^ . 

mai n 


(28) 
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and 


AP 


cool 


32Cf cool L cool( d ") 2 



(29) 


(For the derivation of C, see Appendix D.) 

The pressure drop in the cool down nozzles can be expressed as 



(30) 


When Eq. (27) is solved for the flowrate dm by substituting into it Eq. (28) - 
(30), the following expression is obtained: 


C • f . 
main main 


s RT 
_ + C__ 


meter 


f + c 

cool cool nozzle 


where C. = C . , C , , or C refers to the right hand side of the 

i main’ cool’ nozzle 3 

individual equations without the flowrate term. 

Since the Fanning friction factors in Eq. (28) and (29) also are 
functions of the flowrate, an iterative procedure is used to calculate f . 
and f c The friction factors are initially set to 0.005 and used to 
compute the flowrate. The flowrate is used as the input parameter to FUNCTIO 
FRICT, which computes new values for the friction factors. The new values fo 


f. and f -i are used to recalculate the flowrate, 
main cool 

the friction factors change by less than 0.1 percent. 


This continues until 


After the receiver tank has been chilled to allow two phases of the 


propellant to be present, the nozzle constant, C^, and the length of the cool 
down piping, L coq1 , are set to zero. The contribution of the cooldown line 
and nozzles to Eq. (27) and (31) is negated. This simulates the switching of 


a three-way valve to the closed position in the cooldown line. 
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FUNCTION FRICT 


FUNCTION FRICT calculates the Fanning friction factor corresponding to 
flow through a pipe. Input variables are the flowrate, temperature, density, 
and pipe diameter. 

The Fanning friction factor is related to the Reynolds number and pipe 
diameter. While a variety of functions are available [8], the following 
relationship was chosen: 


f =< 


-4 log 


10 


e/D 

3.7 


5.02 

N Re 


log 


10 


e/D 

3.7 


14.5 

N D 

Re 


(32) 


The pipe roughness e is assumed to be 0.00015 ft (0.004572 cm), the value for 
commercial steel pipe [6]. The Reynolds number is calculated from the equa- 
tion: 


N Re = Du b p/ M 

N _ 4(dni)M 
^Re 7tDp 


(33) 


The viscosity as a function of pressure and temperature for hydrogen 

cannot be calculated with the Fluids Pack thermophysical properties computer 

programs. However, since relatively small changes in hydrogen viscosity do 

not significantly change the friction factor, an average viscosity is used. 

-4 

This value is 1.8 x 10 g/(cm*s) [3] and produces less than 0.5% error in the 
friction factor. 
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FUNCTION PUMPFLO 


FUNCTION PUMPFLO calculates the flowrate between the ET and receiver tank 
and accounts for the use of a pump in the transfer. The overall pressure drop 
can be expressed as 


AP system = AP 11r,e (dn,2) + AP meter + AP cool (dn|2) ' AP pump ( S < 34 > 

See the previous discussion in the section on FUNCTION FLO for descriptions of 

the terms AP . , AP, . , AP . , and AP , . The term AP is a 

system’ line’ meter’ cool pump 

function of the inverse flowrate, 1/dm, and is given by 


AP 


pump 


bR P_H 
dm 


(35) 


The effect of AP „ is opposite that of the other terms in Eq. (34); hence, 
pump 

the opposite sign. 

An initial guess for the flowrate is used to calculate the theoretical 
pressure drop or rise. The flowrate dm is iterated by using a modified Newton 
method [9] which is found to be stable for this application: 


dm n+l = dm n ' 


(AP - AP n )(dm n -dm n _ 1 ) 
< AP „-1 - AP n> 


(36) 


The iteration ceases when calculated and actual pressure drops differ by less 
than 0.01%. 
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FUNCTION PEFF 


The pump efficiency is calculated for an input of the flowrate from FUNC- 
TION PUMPFLO. Two different sizes of pumps are used to transfer the LOX and 
LH 2 propellant. A linear least squares fit of data for the pump efficiency 
versus flowrate [10,11] is used to provide a five-coefficient, fourth order 
equation for each pump. 


The equation for the pump efficiency for LH 2 is: 


H = 0.39213 + 0. 096620(dm) + 2.7793 x 10 _5 (dm) 2 
- 8.1723 x 10" 8 (dm) 3 + 2.5008 x 10" 1:L (dm) 4 . 
The equation for the pump efficiency for LOX is: 

q = 0.32031 + 0. 48261(dm) - 1.6223 x 10 _3 (din) 2 
+ 3.9836 x 10" 6 (dm) 3 - 4.8853 x 10~ 9 (dm) 4 . 


FUNCTION TH 


(37a) 


(37b) 


SUBROUTINE COOL uses FUNCTION TH to calculate a value for the heat of 
expulsion for the vapor phase calculations. This thermodynamic value is de 
fined as: 


6 = \) 



(38a) 


or 



(38b) 


The Fluids Pack subroutines DPDD and DPPT are used to find the partial deriva- 
tives in the numerator and denominator, respectively, and Fluids Pack 
functions are used for the specific heat and density. 
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FUNCTION CWALL 


FUNCTION CWALL returns the heat capacity of the receiver tank wall at a 
particular temperature. It is assumed that the oxygen tank is constructed 
using Inconel X-750* and the hydrogen tank is constructed using aluminum. The 
heat capacity-temperature data for these materials are of the form: 


C 


w 


exp(C 1 + C 2 /T + C 3 /T 2 ) 
1000 


(39) 


The constants C^, and were determined by taking heat capacity- 
temperature data and applying a least-squares technique to the above equation. 
These values are shown in Table 5. The data used are found in Appendix C 
[12,13]. 


Table 5. 

Tank Wall Heat Capacity 

Constants 


Inconel 

Aluminum 

c i 

6.482 

7.293 

C 2 

-106.3 

-121.9 

C 3 

329.1 

347.9 


3 . 5 Program Termination 

SUBROUTINE PFAIL is the normal route for program termination. ENTRY DONE 
is used when the supply tank is empty, i.e., when all the propellant has been 

*Such identification does not imply recommendation or endorsement by the 
National Bureau of Standards. 
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transferred. ENTRY FULL is used when the receiver tank is full. This condi- 


tion occurs if the receiver tank volume is too small, if the quantity of 
propellant is increased with the tank volume remaining the same, or if the 
mass balance on the fluid in tank is wrong. The main entry PFAIL is used to 
indicate a high backpressure condition. This can happen when the receiver 
tank pressure is greater than the supply tank pressure and a pump is not used, 
or when the flowrate is very small. The program will terminate if there is 
improper venting, i.e., when there is a measured flow through the vent, but 
the receiver tank pressure is less than the set vent pressure. 

ENTRY PDATA in SUBROUTINE PFAIL is used when a programmed stop is 
executed. This exit is used in several locations in the main program and cool 
down subroutine, and it is also used by functions DVDPV, DVDPL, DUDPL, DTDP, 
ALPHAL, FLO, and PUMPFLO and subroutine SUPPLY. The COMMON parameter DUMP 
normally has a value of zero until a cause for termination is reached in one 
of the previously listed program sections. DUMP is reset to a value of 1.0. 
When the condition for termination occurs in a function, the program flow is 
returned to either SCAVAG or COOL. The run continues to the end of the present 
integration loop, when the parameter DUMP is evaluated. If DUMP equals 1.0, 
ENTRY PDATA is called, and the output files are replaced prior to termination. 

Several subprograms use convergence criteria to obtain specific thermo- 
dynamic conditions. SUBROUTINE COOL checks for the condition when the density 
in the tank is greater than the saturation density, i.e. both liquid and 
vapor are present in the tank, and the calculations return to the main 
program. If two phases are not obtained within the array size limit, the 
program prints an error message and stops. 

Functions DVDPV, DVDPL, DUDPV, DUDPL, DTDP, and ALPHAL use derivative 
approximations that depend on the change in thermodynamic properties about a 
small interval. The first five functions listed above evaluate the specific 
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volume, internal energy, or temperature about a small pressure interval, while 
ALPHAL evaluates the density about a small temperature interval. In all 
cases, the program will terminate should the convergence of the properties 
about the specified interval not be reached within 100 iterations. The 
message, "KOUNT exceeds 100," noting the particular function, will be printed. 

FUNCTION FLO uses the Fanning friction factor calculated in FUNCTION 
FRICT as its convergence criterion. The flowrate found in FLO is an input for 
FRICT, which in turn determines the friction factor used in the flowrate 
calculation. The convergence of the successive values should occur within 100 
iterations, or the message, "KOUNT exceeds 100 for FLO," is printed on program 
termination. 

FUNCTION PUMPFLO checks the convergence of the actual pressure difference 
between the tanks and the calculated pressure drop through the piping using a 
pump. If the convergence between the two does not happen within 100 itera- 
tions, a message is printed and the program is stopped. 

The supply tank thermodynamics are checked for errors in calculation or 
intolerable operating conditions in SUBROUTINE SUPPLY. When the calculated 
liquid enthalpy is significantly different from the enthalpy found in Fluids 
Pack, the program terminates with the message, "ENTHALPY DISPARITY," and 
prints the two values in disagreement. If the quality of the fluid is greater 
than one part in one million, there exists the possibility for cavitation. 

The program terminates with this condition and an error message, "ET LIQ 
BOILING. STOP TRANSFER." 


3. 6 Programming Information 

A brief description will be given for a typical run of the program. The 

information provided in this section is based on procedures used with the CDC 
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Cyber 750 Computer System and the NOS Version 1 operating system on location 
at NBS in Boulder. The program is written in FORTRAN IV. These procedures 
should be transferable to the computer system at NASA- JSC with minor 
modifications. See Appendix G for program listings. 

File Manipulation 

All computer programs, data, and control procedures for running PROGRAM 
SCAVAG are stored in files. At NBS, SCAVAG is stored in filename COMBO; the 
input data are stored in filenames DATAIN or FUEL; output data are located in 
filenames PSDATA and PSOUT; and the control procedure is stored in filename 
SUBCOM. 

Before running PROGRAM SCAVAG, the operator must prepare the control 
procedure and the input files. The control procedure that NBS has used, 
SUBCOM, retrieves COMBO, DATAIN, and BOXB (containing Fluids Pack) from 
permanent storage and places them into a local or usable state. SUBCOM 
compiles COMBO and, if the compilation is free of errors, loads and executes 
the compiled version. 

The operator enters the following command on a time- sharing terminal to 
run the program: 

SUBMIT, SUBCOM 

This statement is a batch command which means the operator cannot interact, 
alter, or stop the execution of the program or job. To determine when the 
program has stopped, the operator can check the status of the job, or a file 
named DAYF. This file is created by SUBCOM to record the real time sequence 
of control procedures performed by the computer during program execution. 
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Input 


The program operator must determine which propellant and values for the 
initialized parameters are to be used. If LOX is the desired propellant, the 
namelist FUEL must be modified to the following form: 


$FUEL 
NGAS=2 , 

$ 

If LH 2 is chosen, FUEL should have either of the following forms: 

$FUEL 

$ 

or 

$FUEL 
NGAS=1 , 

$ 

Namelist FUEL is contained in a file named FUEL. No decimal point is used 
since NGAS is an integer. 

To change values for the initialized variables, similar steps are used 
with namelist NAME. NAME is contained in a file called DATAIN. Any number of 
changes can be made for the parameters listed in Table 1. For instance, to 
change the initial supply tank temperature and pressure, the quantity of 
propellant to be transferred, the pump power, and the length of straight pipe 
for the LH 2 system, one would alter NAME to the following form: 


$NAME 

TET = -400. , 

PET = 25. , 

TOTAL = 5000. , 

POWER = 2. , 

LENGTH = 100. , 

$ 

Since all the chosen variables are real, the values must have decimal points. 
All values are entered in English units. If no changes are desired, NAME must 
have the following form: 


$NAME 

$ 



The program statement defines TAPE5 and TAPE7 for use during the program. 
SUBCOM, which directs the sequence of operating steps for running the program, 
sets TAPE5 equal to filename DATAIN, or namelist NAME, and TAPE7 equal to 
filename FUEL, or namelist FUEL. These files must exist in order to run the 
program. 

Output 

The output from program SCAVAG is directed to two locations. While SCAVAG 
is running, the graphic display data is written to TAPE8 and stored in filename 
PSDATA. When SCAVAG terminates, the termination message and initialized par- 
ameters are written to TAPE6 and stored in filename PSOUT. See Table 6 for 
sample output for LH 2 . The operator can access PSOUT at a later point for 
further evaluation. 

Two FORTRAN programs utilize the data in PSDATA to produce plots des- 
cribing the transfer. The first program is called PLOT. At NBS PLOT is 
attached to the job by the control procedure SUBCOM and run after SCAVAG has 
terminated. The plots that are drawn use the computer systems mathematics 
library (called STARPAC at NBS in Boulder) and are written to TAPE6, following 
the output described above. The entire contents of TAPE6 are stored in PSOUT. 
(Further explanation is given in Appendix F.) 

The second program which uses PSDATA is called DEMO. The graphics pack- 
age called DISSPLA (15) is utilized to produce plots which are more refined 
than those generated through PLOT. (Additional details about DEMO can be 
found in Appendix F.) The plots are directed to one of three output devices 
by selecting the appropriate CALL statement in DEMO. For the computer system 
at NBS CALL HP7221 accesses a multiple pen plotter; CALL TK4010(960) accesses 
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Table 6. Sample Output for Hydrogen 


Supply tank initial conditions: 


Pressure 

32.000 

psia 

(0.221 MPa) 

Temperature 

-425.000 

F 

(19.1 K) 

Hydrogen Mass 

3098.00 

Pounds 

(1405.2 kg) 

Heat Leak 

90000.000 

Btu/hr 

(26376.3 W) 

ver tank initial conditions: 

Pressure 

18.000 

psia 

(0.124 MPa) 

Temperature 

-290.000 

F 

(94.1 K) 

Vent Pressure 

30.000 

psia 

(0.207 MPa) 

Tank Volume 

780.000 

Cubic Feet 

(22086.5 kg) 

Tank Wal 1 Mass 

975.000 

Pounds 

(442.2 kg) 

Heat Leak 

0.000 

Btu/hr 

(0.0 W) 

fer Parameters: 

Pipe Diameter 

5.000 

Inches 

(12.7 cm) 

Length of Straight Pipe 

70.000 

Feet 

(21.34 m) 


10 Elbow(s) 

0 Gate Valve(s) 

2 Globe Valve(s) 

0 Angle Valve(s) 

0 Butterfly Valve(s) 

1 Flow Meter (s) 

Heat Leak Into Piping 
Pump Power 


50400.000 Btu/hr 
0.0 HP 


(14770.7 W) 
(0.0 W) 


Cool-Down Parameters: 


Length of Piping 
Number of Nozzles 
Nozzle Diameter 
Header Diameter 


20.000 Feet 
6 

0.500 Inches 
2.000 Inches 


(6.10 m) 

(1. 3xl0- 2 m) 
(5. Ixl0- 2 m) 


Cool-Down Time 


0.1754 minutes 


At 2.344 minutes from start of propellant transfer the supply tank is empty. 

At this time the, receiver tank is at -422.89 °F (20.3 K) , 16.23 psia (0.112 
MPa), and contains 3098.00 pounds (1405.2 kg). 

During the simulation 0.00 pounds (0.0 kg) were vented. 

Supply tank is at -424.9907 °F (19.1 K), 31.55 psia (0.218 MPa). 
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a graphics terminal; CALL FR80(0.0) accesses microfilm. The running of 
program DEMO is separate from running SCAVAG and PLOT. Examples of plots from 
program DEMO are found in Appendix G. 
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4. RESULTS AND DISCUSSION 


The feasibility of low-g propellant transfer is analyzed parametrically 
with the computer model. Particular emphasis is placed on the thermodynamic 
changes to the system caused by varying the temperature, pressure, or heat 
leak at any point. The principal parameter in this investigation is the 
overall time of transfer, since the transfer is to take place during a mated 
coast period after MECO. Other parameters of interest in the design of such a 
system are the pump size and piping requirements. 

Nozzles are used during the cooldown of the receiver tank to induce a 
flowrate of propellant lower than the flowrate during the main part of the 
transfer. The sizes of the associated piping and nozzles, which are simply 
small orifices situated around the header, have been parametrically determined 
to accomplish a low flowrate without causing an excessive pressure 
drop. This specification holds except in those cases where the transition 
from the main piping to the cooldown piping is great. (Further explanation is 
provided later.) When two phases occur in the tank, the nozzles are shut off; 
the transfer via the main line resumes for the remainder of the simulation. 

Because the transfer can require a pump, the liquid from the ET must be 
subcooled to prevent cavitation. This is most notably the case for LOX which 
has insufficient ullage pressure to accomplish a complete pressurized 
transfer. Single-phase flow from the supply tank also requires sufficient 
thrust to settle the liquid against the outlet port. We assume such a thrust 
for the purposes of this modelling effort. 

When the program is run with the default parameters, the time required to 
complete the transfer of either propellant is relatively short. Transfer of 
LOX with a moderately sized booster pump takes 1.94 minutes; pressurized 
transfer of LH 2 requires 2.34 minutes. If the MECO mated coast period is 
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scheduled to continue to 20 minutes, the scavenging of both cyrogenic 
propellants can be completed in less than 12 percent of the alloted time. The final 
system conditions are shown in Table 7. Plots of the transfer are found in 
Appendix G. 


Table 7. 

Final Conditions 

Using Default 

Values 



0 2 


h 2 


Time of Transfer, min 

1.936 


2.344 


Receiver Tank 

Temperature, °F (K) 

-312.35 

(81.69) 

-422.89 

(20.28) 

Pressure, psia (MPa) 

6.72 

(0.046) 

16.23 

(0.112) 

Mass in Tank, lb (kg) 

6270. 

(2846.6) 

3098. 

(1405.2) 

Mass Vented, lb (kg) 

0 . 


0 . 


Supply Tank 

Temperature, °F (K) 

-314.9998 

(80.22) 

-424.9907 

(19.12) 

Pressure, psia (MPa) 

19.88 

(0.137) 

31.55 

(0.217) 


4. 1 Heat Leaks to the System 

The conditions in the ET are governed by two opposing influences: outflow 
and heat leak. In the absence of heat leak, outflow gives an adiabatic expan- 
sion of the vapor and a resulting decay in temperature. In the absence of 
outflow, heat leak gives an increase in temperature and pressure. Figures 4 
and 5 show the expected results of varying the amount of heat to the ET. 

Figure 4 shows that a large rate of heat flow into the LOX ET must be present 
for a significant change in the system parameters to occur. 


% 
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Figure 4a. 
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Figure 4b. Effect of varying the heat leak to the ET, LFL (default: 
90,000 Btu/hr (2.6xl0 4 W)) 
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Figure 5a. Effect of varying the heat leak to the transfer piping, LOX 
(defualt: 104,400 Btu/hr (3.0xl0 4 W» 
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Figure 5b. Effect of varying the heat leak to the transfer piping, LH ? 

(default: 50,400 Btu/hr (1.5xl0 4 W)) ^ 
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6 5 

The value at which this occurs is on the order of 10 Btu/hr (2.9 x 10 W). 

The pressures in both tanks undergo rapid increases, while the temperatures 
show slight increases. 

The changes in temperature and pressure that the LH 2 system experiences 
are less abrupt than those of the LOX system. They do show the same trends, 
however. As the heat leak to the ET is increased, the pressures and 
temperatures show gradual increases. 

As the heat leak to the transfer piping is increased, the final 
conditions in the receiver tanks for both propellants show a corresponding 
rise greater than the results for a heat leak to the ET. The results are 
presented in Figure 5. Venting occurs for LOX when the heat input is 
increased to 5 x 10 Btu/hr (1.5 x 10° W); approximately 28% of the propellant 
is vented. The LOX temperature of -284.29 °F (97.28 K) at the vent pressure 
is located at the liquid-vapor boundary. The pressure history for this 
transfer is shown in Figure 6. 

Similar observations can be stated for the case of LH 2 . When the heat 

S 5 

leak to the transfer piping is increased to 10° Btu/hr (2.9 x 10 W), a small 
amount of fuel is vented. This occurs at the end of the transfer, as the 
fluid reaches the liquid-vapor boundary. 

Conditions in the ET do not vary with an increased heat leak to the 
piping. This is due to the flow of cryogen from the tank which directs any 
potential effects to the receiver tank. The overall time of transfer 
increases linearly for both LOX and LH 2 . The time required to transfer LOX 
remains quite brief, but the transfer requires an additional minute when 
5 x 10^ Btu/hr (1.5 x 10^ W) leaks into the LH 2 piping. 

The effect of heat leak into the receiver tank is similar to that into 

the transfer piping. For LH 2 the transfer time is increased by nearly a 

• 5 5 

minute for a heat rate of 4 x 10 Btu/hr (1.2 x 10 W). Increasing the heat 
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Figure 6. Pressure history of LOX transfer for a heat leak of 2x10® Btu/hr 
(5.8x10 s W) to the transfer piping 
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6 5 

leak to the receiver tank at the rate of 10 Btu/hr (2.9 x 10 W) causes LH ? 
to vent. The transfer time for LOX rises smoothly with increased heat input; 
an additional minute is required to transfer LOX at a heat leak of 1.2 x 10 ^ 

5 

Btu/hr (3.5 x 10 W). Also, venting takes place at this rate of heat 

addition. These results are shown in Figure 7. 


4.2 Receiver Tank Temperature 

Since the temperature for the LOX receiver tank is unspecified in the 
study by Brux and Stefan [1], the effects of several temperatures were 
examined in this study. The results are shown in Figure 8. It can be assumed 
that the receiver tank will receive some soakback heat leak from the 
insulation as well as radiation prior to MECO, Thus, temperatures higher than 
cryogenic are realistic and allow examination of the chi 11 down process. A 
wide range of initial temperatures neither significantly alters the time 
required to transfer the oxygen nor the final receiver tank temperature. The 
coolest temperature input of -100 °F (199.7 K) would require the vessel to be 
prechilled on the ground. For temperatures above 120° F (321.8 K) a small 
amount of venting (less than 1 percent) takes place at the beginning of transfer. 
This rapid increase of the RT pressure to the vent pressure is shown in Figure 
9. 

The study by Brux and Stefan [1] recommends that the LH 2 receiver tank be 
prechi lied to -290 °F (94.1 K) to prevent venting during the transfer after 
MECO. They state that at higher temperatures, a large amount of venting would 
be required to achieve chi 11 down of the receiver tank. We find that only 
minor venting at the start of the transfer takes place when given an initial 
temperature between -250° F (116.3 K) and -50 °F (227.4 K). For all cases, 

the overall time of transfer is less than 2.80 minutes. 
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Figure 7a. Effect of varying the heat leak to the receiver tank, LOX 
(default: 0 Btu/hr (0 W)) 
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Figure 7b. Effect of varying the heat leak to the receiver tank, LH^, 
(default: 0 Btu/hr (0 W)) 
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Figure 8a. Effect of varying the initial receiver tank temperature, LOX 
(default: 60°F (288.6 K)) 
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Figure 8b. Effect of varying the initial receiver tank temperature, LH„ 
(default: -290°F (94.1 «)) 1 
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Figure 9. Pressure history of LOX transfer for an initial receiver tank 
temperature of 200 °F (366.3 K) 
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4.3 ET Temperature 


Since the previous study [1] indicates the ET fluid is subcooled, the 
default initial temperatures in this study were between the melting line and 
the liquid-vapor boundary [3,4,5]. Increasing the ET initial temperatures 
cause longer transfer times, while additional cooling accomplishes the 
transfer faster. As the melting line temperature of either propellant is 
approached, failure of the Fluids Pack [2] routine occurs. When the 
temperature reaches the vapor-liquid boundary, the simulation stops due to 
vapor induced cavitation. Because the transfer of LH 2 is pressure induced, 
the transfer times are significantly affected by warmer ET temperatures. 
Additional cooling time is required before the switch to the transfer mode 
through the main piping. 


4. 4 Pipe Size and Pump Requirement 

When the values for line size and pump power suggested by Brux and Stefan 
[1] are used in the model, the time required to transfer the nominal quantity 
of either propellant is less than 3 minutes. Because the mated coast period 
after MECO probably will be somewhat longer than this, transfer times of 
approximately 8 minutes or less appear reasonable. Such transfer times are 
obtained with a number of cases where the transfer pipe diameter or the pump 
size is changed. Figure 10 illustrates these changes. 

If the LOX default pump size remains unchanged at 4 HP (2982.8 W) while 
the line size is reduced to 2 inches (5.1 cm), the transfer time is slightly 
less than 7 minutes. When the pump size is reduced by half, the transfer time 
is less than 4 minutes for a 3- inch (7.6 cm) diameter pipe. Line sizes of 
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approximately 1 inch (2.5 cm) appear to restrict the flow too much and cause 
venting in the receiver tank and extremely long transfer times. A 6-inch 
(15.2 cm) line size promotes a fast transfer of the remaining LOX propellant, 
but it also is heavy; no further study of this line size was done. Since the 
smaller line sizes would be lighter, easier to design around the existing 
equipment, and less costly, reducing the line size from the default value of 4 
inches (10.1 cm) to 3 inches (7.6 cm) is a reasonable exchange for the in- 
crease in transfer time. A pump size of 2 HP (14914.4 W) would not tax the 
energy requirements in the system. 

Similar results are seen for LH 2 scavenging. A reasonable transfer time 
of slightly more than 7 minutes is obtained when a 3- inch (7.6 cm) line size 
is used without pump assist. However, there appears to be too much flow re- 
sistance in smaller line sizes, even with the assistance of a pump. Hence, 
longer transfer times are seen for these cases, even when a large pump is 
used. When a 2 HP (14914.4 W) pump is used with 3-inch (7.6 cm) piping, the 
time is reduced to slightly less than 6.5 minutes. Reducing the line size to 
3 or 4 inches (7.6 or 10.1 cm) and pressurized transfer of LH 2 is an optimum 
combination for the system. 


4. 5 Limitations 

Accuracy for the results is limited by the size of the time step used in 
program. When a very small step is used, the amount of computer time required 
is large and consequently quite costly. Thus, the time steps used in the 
previously discussed analyses are large enough to allow speedy computation, 
yet small enough to yield reasonable accuracy. When the time step is reduced 
by a factor of seven in SUBROUTINE COOL and by a factor of thirteen 
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in SCAVAG, the total time required to transfer LOX propellant changes by less 
than 1.5 percent. However, the computer time required for this run increases by a 
factor of 80. 


5. SUMMARY 


We have developed a computer model to investigate the scavenging of cryo- 
genic propellant from the ET to shuttle mounted receiver tanks in the cargo 
bay. Based on the requirement that sufficient thrust exists to settle the 
liquid in the ET, the model predicts the successful transfer of liquid oxygen 
and hydrogen in short periods of time, typically less than 2.5 minutes for 
both LOX and LH 2 using the default parameters in the program. A parametric 

analysis was performed to evaluate the effects of heat fluxes to various 

system components, prechilling the receiver tank, and line and pump size 
variations. 

Moderate to large heat fluxes affect the pressurized transfer of LH 2 more 

than the pump-assisted transfer of oxygen. Heat leaks to the ET increase the 

pressure, and thus speed up the transfer of LH 2 . Heat leaks to the piping or 
receiver tank slow down the transfer of both propellants due to additional 
receiver tank chilling demands, and cause venting at sufficiently high fluxes. 

Prechilling the receiver tank for LH 2 is a critical requirement. The 
LH 2 receiver tank must be approximately -290 °F (94.1 K) in order to transfer 
the fluid without venting. The LOX receiver tank can be close to 300 °F 
(421.9 K) without the significant venting of fluid. 

Line and pump size selection are dependent on the time allowed, the 
piping design constraints and energy demands of the transfer system. 
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Pump-assisted transfer in the LOX scavenging is needed, and by using a 2 HP 
(1491.4 W) pump and a 3-inch pipe diameter, a relatively lightweight, fast 
system can be used. Liquid hydrogen can be transferred by pressure difference 
in a 3-inch line. The smaller line and pump sizes allow the transfer to be 
completed in less than 8 minutes. 
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APPENDIX A 


Derivation of single-phase receiver tank equations: 

First Law 

dU = '“’total + h l d ”l * h 2 dm 2 
Mass Balance dm = dm^ + dm 2 

By applying the mass balance to the system and assuming h = h 2 , the first law 
may be expressed as: 

mdu + udm = dQ totg ^ + h^dm^ - h(dm^ - dm) 

. . '“’total . h l dm l h<d V dnl) udm 

du - + - - 

mm mm 

du = + ir < h i - w + nr 2 • (a.i; 


Because we wish to find temperature as a function of mass flow we next 
expand du. The total differential energy equation with respect to temperature 
and density is given by 


du = ^ dT + p 


9u\ dg 
9p/ j P 


du - C dT - v f !rA ^ • 


\> ' ' \3v/y p 


(A. 2) 


By using the thermodynamic relationship 
du = Tds - Pdv 
we can obtain , . 

* (i) T ■ * (§4 - - ■ 


The definitions for and the Gruneisen parameter <t> are 


C v = T l 


3s) 

si) 


v 


fifel 


T v®P 
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(A. 3) 


and the chain rule gives 

“(fi = * TC 


v P 


Equation (2) then becomes 


du = C dT - (d)TC 
v Y v 


P-. 

p' m 


Combining Equations (1) and (3) and solving for dT yields: 

C dT - OTC - -) — = (h - h) + 

v VY vp / m m m v l y m 


or 


dT = - 
m 


„ TJ . < h l - h \ , . d ^total 

(tj)T dm + p ) drn-. + ~ 


(A. 4) 


v v 

The heat flux is the sum of the heat leak through the tank and the heat 
in the wall of the tank, 

dQ total = + d(! w 


dQ = - m C dT 
w w p w 

Substituting these expressions into the temperature equation yields the 
following expression: 

dT ♦ "« C Pw dT = i f |T| „, . < h l- h > dm. ♦ ' 

< r-—~ m UTdm + ~ 1 r } 

C v m V C v V 


Upon rearrangement the following relationship is obtained: 


dT = 


()>Tdm + ( h-^ -h)dm^/C y + dQ/C y 

m + m w c p„ /c » 
w 


(A. 5,3) 


We next wish to find pressure as a function of mass flow. The total 
differential energy equation with respect to pressure and density is given by 

du = (1) p (|u) dP ♦ p (!) p * . 
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The first term on the right-hand side contains the Griineisen parameter 


which is also defined as 


• ■ m p ■ 

The second term on the right-hand side may be expressed with the 
expulsion 


6 = - P 


ah\ 

3p/r 


and by using the definition of enthalpy to give 


/3u\ /3h\ P 

P \3p/p P \3p/p P 


The differential energy equation may then be written as 

du = \ dP - (0 - -) ^ . 

P4> V P / P 

Combining Equations (A.l) and (A. 6) and solving for dP when dm^ 
yields the following expression: 


dP = i ^ [0 + h 1 - h] dm + 4>(dQ - m w C w dT)j . 

When dm^ is not equal to dm, there is flow through the vent. Thus, dP 


dm 


[(h - h 1 )dm 1 - dQ + m^dT] 


heat of 
(A. 6) 


(A. 7) 
dm 

(A. 8) 
= 0 and 

(A. 9,5) 
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Appendix B 

Derivation of two-phase receiver tank equations: 

First Law 

dl) = dQ - dW + I h-jdm-^ 

Since there is no work done, the first law may be expressed as: 

'A + m v du v "• u J- dn V + u v dm v = ‘“’total + h l d, "l + h 2 dm 2 (Ba) 


Continuity 


dm^ + dm v = dm^ + dm2 


(B. 2,7) 


Constant volume 


V m £ + + V v dm v + m v dv v = 0 • 


(B. 3,9) 


Solving for dm 


& 


v , 

dm t = dl \ - ^ + "v du v> 


Let 




The equation for the differential mass of the liquid in the tank is: 

v. 


d m = - — ~ dm - — 
£ V 


m 1 9V \ 4- m i 9V 

v Up/ v j 


dP 


(B. 4) 


Substituting Equation (B.4) into Equation (B.2) produces the following expres- 
sion: 


y 

— - dm + dm - — 
V £ V V 


9v \ l 9v\ 

"z V 9 p/ £ M^/v 


dP = dm-, + dm,. 

1 t 
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Solving for the differential mass of the vapor: 


dm. 


dm^ + dm 2 + 


'l 


‘ 3v ' 


3v 


\ 9 P A + m vWv 


(1 vV 


dP 


(B. 5) 


Substituting Equation (B.5) into Equation (B.4), we obtain the following 
equation for dm £ 


dm. 


(v v /v £ )(dmi + dm 2 ) - v £ 


3v>\ + /9v) 

"i l 3p A m vV8( Vv 


dP 


(i - * v /V 


(B. 6,10) 


Combining Equations (B.l), (B.5) and (B.6), we obtain 


„ , V rtp + „ (*>\ HP „ -!y “"l 4 dm 2 

MapJ/P ”v\3pj/ P U £ v t (1 - » v /V 


§v\ / 3v\ 

V 3 P/j + "»V 8 7v df \. 

+ (1 - VV — 


dm^ + dm 2 
u v (1 - v v /v £ ) 


1 / 9v | / 3\> 

v £ m A 9p k + % W v dP 

(i - VV 


dQ total + h l dm l + h 2 dm 2 


(B. 7) 


We can solve Equation (B.7) for the differential pressure when there is no 
venting: 


dP 


h l “ u v + ^ 


V 


"fTpT 


dm l + d( Wal (1 " W 


. (B. 8,11a) 
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where 


F(p) = m £ 



Ujt ) 



- m C 


w w 


y 3P ysat 




(B. 9,11b) 


When there is venting, the pressure is constant. Since dP = 0, the following 
equation is obtained when Equation (B.7) is solved: 

( v dm^ + dm2 \ / dm^ + dm^ \ 

(1 “ W/ \ Uv (1 - Wj 

+ dQ - m C dT + t^dm., + h 0 dm 0 . 
w w i i c c 

Finally, an expression for the differential mass through the vent is obtained: 


£ v 


dm. 


£ 


U v + hj (1 


V, 


V, 


dnt-j + dQ (1 


£ 


(h, 


h £^ 


(B. 10,12) 
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SPECIFIC HEAT. J/(kg K| 


Appendix D 


Conversion factor calculation 

The conversion factor C for Equation (26) may be derived as follows 


where 


Therefore, 
32C = 


„ r . AP p7i 2 D 5 

u fL(dm) 2 M 


AP 

P 

D 


atm 

gmol/£ 

cm 


L = m 
M = g/gmol 


1 kg 


1 SL atm 


lm 


Y 


(atm)(gmol/£)(cm 5 ) 

(m)(gmol/s) 2 (g/gmol ) ^100 cm^ 

. ( 1000 gVlQl.325 j \ 


1 SL 


0.001 m 3 


1 kq m2 
J s 2 


32C = 3.158 


(atm)(s 2 )(cm 5 ) 

(m)U)(g) 
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Appendix E 


Derivation of the supply tank equations: 


First Law 


dll = dQ - dW + Zh.m. 

The first law can be written for both phases, and when the mass of vapor is 
held constant we have the following expression: 


vapor phase m v du v = d Q v " p ( m v dv, v ) 


Since 


m v h v = Vv + P Vv 


dh - v dP = du + Pd-u . 
v v v v 


Upon substitution and rearrangement, the vapor phase enthalpy equation becomes: 


dh = — i + v dP 
v m v v 


(E. 1,13a) 


For the liquid the first law may be expressed as 


liquid phase m^du^ + u £ dm £ = d Qj? ~ Pm £ dv £ “ pv £ dm y + ^£ dm A 


By using the previous substitution and rearrangement, the liquid phase enthalpy 
equation can be written as 


dh i = 5J- + V p 


(E. 2,13b) 


Since the £7 volume is constant the following expression is used: 


V = constant = v 0 m 0 + \> m 

X, Xs V V 
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Upon rearrangement. 


v^d'V = " ^ m £ dv £ + m v dv v^ ' ( E -3) 

The total derivative for volume with respect to pressure and enthalpy is 
expressed as 

dv = (W) p dh + (sp) h dP ' *- E ' 4 ^ 

We recall the definition for the heat of expulsion 0 and substitute it for the 
first term in the total volume derivative as follows: 



To elucidate the thermodynamic nature of the second partial derivative, 
the following thermodynamic identity is used: 

S), (8), (8), ■ - ■ ■ 

The term^ |jj^ can be expressed in terms of 8 

3v\ / 9v \ ( 3v\ _ v / 9h\ 

3P/ h \ 3h /p \3P/ V ~ 8 V 3p / v ‘ 


The definition for the Grilneisen parameter can be used with the first law to 
eliminate the remaining partial derivative. The Griineisen parameter is: 



By writing the first law in terms of derivatives with respect to pressure at 
constant volume, we can make the following substitutions: 
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9h\ J 9u 


ap; v ^ap/ v 

= ^ + V 
<!> 


( 1 ) + 1 


<e r 


The second term in Equation (E.4) can be expressed as: 


v /v(<t> + 
8 ( (J) 


Equation (E.4) can now be expressed by the following relationship 


dv = l dh - At* 1} dP . 

0 0ffl 


We can substitute Equation (E.5) into (E.3) to obtain the following 


relationship: 


A ^ Ii ( j + i \ 

\ e * V« s . e i ' ♦ k 


VA 


By solving Equations (E.l), (E.2), and (E.3) simultaneously, a relation- 
ship for the ET pressure is found: 

dp . »|V, (£ 7 14 . 

, V? ' 
e A VtVv 


An expression for the temperature of the liquid can be obtained by using 


the liquid enthalpy equation (E.2) and the total enthalpy derivative. 
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The total enthalpy derivative with respect to temperature and pressure is 


given by: 


dh 



dP + 



(E.8) 


The second partial derivative in Equation (E.8) is the definition of the con- 
stant pressure heat capacity C . The first partial can be expressed as [14] 

P 



= v 



The definition for the bulk thermal expansivity a of a fluid is given by [7] 


a = 


If §£\ = I ( §v 

p\8T/ p v \3T 


The partial derivative may be expressed as 


(§f) t = ” (1 ■ °) • 

By substituting for the partial derivatives as discussed, Equation (E.8) may 
be written as follows for the ET liquid phase 


dh £ = v (1 - or A ) dP + C p dT . 

By substituting Equation (E.2) for dh £ in the previous equation and 
solving for dT, the temperature equation is given by the relationship 


1 dQ £ 

dT = r + V£ dP • (E - 9 ’ 15) 

p ^ 
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Appendix F 
Plotting routines 


The data which describe the transfer of fluid are stored in file PSDATA 
(see printout of file in Appendix C). These data are used in two plotting 
programs called PROGRAM PLOT and PROGRAM DEMO. Both programs are written in 
FORTRAN IV and are fully described below. The data file is described first. 

FILE PSDATA/TAPE8 

The first four lines of the file are data points for identifying and 
initializing the plots for either PROGRAM PLOT or DEMO. The first line of 
data contains four identifying parameters -- the transfer line size, pump 
size, and the initial receiver tank temperature and pressure. The second line 
contains the alphanumeric name of the propellant, HYDROGEN or OXYGEN. The 
total number of data lines, N, is written on the third line. Finally, the 
fourth line holds the lower and upper limits for each of the plots to be 
drawn. The order of these pairs of values is the time, receiver tank 
temperature and pressure, the flowrate, the mass of vapor, mass of liquid, 
total mass of fluid transferred, the quality of fluid in the ET, and the 
quality of fluid entering the receiver tank. 

Following these lines are the N lines of data. Each entry holds the data 
describing the transfer for each increment of time. The values in each line 
are the elapsed time, the receiver tank temperature, receiver tank pressure, 
the flowrate, the mass of vapor, liquid, and total mass of propellant 
transferred, the ET fluid quality, and the quality of fluid entering the 
receiver tank. 
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PROGRAM PLOT 


The NBS computer system in Boulder has a mathematics library called 
STARPAC from which subroutines may be called to perform assigned tasks. For 
program PLOT the SUBROUTINE PPC is used to produce four plots and is accessed 
by the statement CALL PPC(A1, A2, A3, A4, A5, A6, A7, A8, A9, A10). Each 
parameter is described as follows: 

A1 Array containing the value to be plotted on the vertical axis. 

A2 Array containing the values for the elapsed time during the 

transfer to be plotted on the horizontal axis. 

A3 The number of lines of data to be read. 

A4,A5,A6 Parameters describing whether the plot is to be logarithmic, 
the size of the plot is to be altered, or the plot is to be 
reprinted. All are set equal to zero. 

A7,A8 Lower and upper values for the vertical axis. 

A9,A10 Lower and upper values for the horizontal (time) axis. 

PLOT first assigns single array dimension sizes for the time, tempera- 
ture, pressure, flowrate, mass of vapor, and mass of liquid. Second, the 
first four lines of data are read into appropriate variable names. Following 
these assignments the data points are read into the dimension arrays. 

The first call to PPC enables a plot to be drawn of the receiver tank 
temperature as a function of time. The second call to PPC is for a plot of 
the receiver tank pressure as a function of time. The third call plots the 
mass of propellant in the vapor phase as a function of time. The fourth call 
to PPC produces a plot of the mass of liquid phase propellant as a function of 
time. Each plot has an appropriate heading. There are 10 vertical and hori- 
zontal axis divisions for each plot. 

When batch control procedure SUBCOM is used, the plots are located after 
the output from program SCAVAG and all output subsequently are written to 
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TAPE6 and PSOUT. Should the operator choose to have the output separated, two 
other batch control procedures are used. First, SCAVAG is run with NOPLOT, in 
the same manner as has been described in the main text for SUBCOM. Output is 
written to TAPE6/PS0UT and PSDATA as before. To generate the plots for the 
transfer with program PLOT, control procedure LIBPLOT is used. Output is 
directed to TAPE6 and stored in filename TRACE. (See Appendix G for the pro- 
gram, data, control procedure listings, and sample plots.) 

PROGRAM DEMO 

The DISSPLA graphics package [15] is available at NBS in Boulder and at 
NASA- JSC. Since it produces high quality plots and is quite versatile, 

DISSPLA was chosen for graphic display of some of the results of this study. 

It is used by program DEMO to produce seven plots. Two plots, each allowed 
half of a standard page, are placed on the first three pages, and the re- 
maining plot is located on the upper half of the fourth page. 

DEMO assigns dimension sizes and variable names in a manner very similar 
to that of program PLOT. However, additional arrays and variable names are 
required for the parameters used in the DISSPLA subroutine calls. After the 
lower and upper bounds for each of the plots are read from PSDATA, DEMO calcu- 
lates the division size for the axes of each plot. After the data points are 
read into their respective arrays, the program is fully initiated and proceeds 
with a succession of DISSPLA subroutine calls which set up and plot the curves. 
A full description of each subroutine can be found in the DISSPLA user's 
manual. (The program listing and sample plots are located in Appendix G. ) 

Several items are unique to program DEMO. Whereas PLOT reads the first 
line of data into four dummy variables, DEMO uses these data as identifiers 
located in the corner of the first plot. The name of the propellant is not 
used in PLOT, but in DEMO it is used in the heading for each page. In 
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addition to the four plots which are generated by both PLOT and DEMO, DEMO 
also plots the total mass of propellant transferred, the flowrate into the 
receiver tank, and the percent of vapor and liquid phases in the receiver 
tank, as a function of time. 

DEMO is run interactively by the operator. Instead of storing the output 

in files which are accessed at a later time (i.e., batch operation, SCAVAG 

with control procedures SUBCOM or NOPLOT, and PLOT with control procedures 

SUBCOM or LIBPLOT), output is generated immediately at the terminal or other 

output device. To run DEMO at NBS, the operator types in at the terminal 

GET, DEMO 
DEMO 

or BEGIN, DEMO, DEMO 

The first line instructs the computer system to retrieve program DEMO from 
storage and place it in the local mode and ready for use. When DEMO is typed, 
the procedure for running the program is activated. The alternative command 
identically performs both statements. The procedure is located in the lines 
prior to the program statement and uses a sequence of control statements 
similar to those contained in the batch control procedures SUBCOM, NOPLOT, or 
LIBPLOT. 

Output is generated at the terminal screen when the appropriate call is 
made and graphics capabilities are available. Other output options are a 
multiple pen plotter or microfilm. The first page, or two plots, are drawn 
automatically. The operator types 

COPYBR, Z 

to have each subsequent page drawn at the terminal or the pen plotter. All 
plots are automatically drawn on microfilm. 
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Appendix G 

Program listings and sample plots 
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PROPELLANT SCAVENGING PROGRAM-COMBO 


PROGRAM $CAVAG(INPUT#TAP£5* TA PE 6# TAPE7* OUTPUT »TAPE6*TAPE8) 

THIS PROGRAM SIMULATES THE FILLING OF A RECEIVER TANK IN LOW G 
USING A SUPPLY/EXTERNAL TANK. 

COMMON/ PAR AM/ MW# MW ALL # 0/ VOL> OP I PE# OET * TOT AL # P VENT 
COMMON/ SUB/ FLOWl* FL0W2#PR£SS*TEMP#TIME*HET 
COMMON/FLOW/ DIAM,0ET#TET*P£T*LEN6TH,PMETER 
COMMON/NOZZLE/CDIAM#NOZDIA#H£ADIA»NN0Z*CONST#CDLENG 
COMMON/PUMP/ POWER* EFF 

COMMON/ PLOT/ NPL0T,X*Y1,Y2#Y3*Y4*Y5# Y6*Y7*Y8 
COMMON/PROP/NGAS 

COMMON/ETANK/ETLIO*ETVAP*VET»ETO 
COMMON/ PROB/DUMP 

REAL MVAP,HLIO,MV*MVALL*L£NGTH*NOZDIA 

REAL X(600)#Y1(600)*Y2(600)#Y3(600)*YA(600)#Y5(600)#Y6(600) 

♦ # Y7(600)#Y8(600) 

THE NAMELIST "NAME" CONTAINS ALL PAR AMFTERS THAT CAN BE VARIED 
FROM DEFAULT VALUES. 

NAMEL I ST/NAME/ PET# TET* OET* TOTAL *DI AM, NE LBOW* NANGL E# NG ATE# 

♦ NGL0BE#NBUTT#LENGTH,METER*QPIPE*P0WER# VOL *MWALL# TEMP* PRESS# 

♦ PVENT,Q,CDLENG»HEA0IA#NN0Z*N0ZDIA#VET»C0NST 

NG AS « 1 FOR HYDROGEN AND 2 FOR OXYGEN CHANGE IN NAMELIST "FUEL" 
NAMFLIST/FUEL/NGAS 

NG AS ■ I 

RE ADI 7* FUEL ) 


INITIALIZE FLUID PROPERTIES WITH FLUIDS PACK SUBROUTINE. 
IF(NGAS.E0.2)G0 TO 1 
CALL DATA P H2 

INITIALIZE MOLECULAR WEIGHT 

MW » 2.0159* 

SET DEFAULT VALUES FOR THE SIMULATION. 

INPUT PARAMETERS FOR SUP B LY TANK t 
PET • TANK PRESSURE# PSIA 
TET - TANK TEMPERATURE# DEG. F 
OET ■ HEAT FLUX INTO SUPPLY TANK* BTU/HR 
TOTAL » TOTAL PROPELLANT IN SUPPLY TANK# POUNDS 
VET • TANK VOLUME* CUBIC FEET 


PET 

m 

32. 

TET 

m 

-425. 

GET 

m 

90000. 

TOTAL 

- 3098 

VET 

m 

53518. 


INPUT PARAMETERS FOR PIPING BETWEEN TANKS* 
DIAM « PIPE DIAMETER# INCHES 
NELBOW - NUMBER OF RIGHT-ANGLE ELBOWS 
NGATE - NUMBER OF GATE VALVES 
NGLOBE ■ NUMBER OF GLOBE VALVES 
NANGLE - NUMBER OF ANGLE VALVES 
NBUTT « NUMBER OF BUTTERFLY VALVES 
LENGTH - LENGTH OF STRAIGHT PIPE# FEET 
METER • NUMBER OF FLOW METERS 
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OPIPE « HEAT FLUX INTO PIPE, BtU/HR 
POWER - PUMP POWER, HP 

DIAM ■ 5. 

NELBQW - 10 
NG ATE • 0 
NGLOBE • 2 
NANGLE « 0 
NBUTT • 0 
LENGTH > TO. 

METER ■ 1 
OPIPE • 50400 • 

POWER • 0. 

INPUT PARAMETERS FOR RECEIVER TANK * 

VOL ■ TANK VOLUME, CUBIC FEET 
MWALL • TANK WALL MASS, POUNDS 

TEMP « INITIAL RECEIVER TANK TEMPERATURE, OEG F 
PRESS • INITIAL RECEIVER PRESSURE, P$IA 
PVENT - RELIEF VALVE VENT PRESSURE, PSIA 
0 » HEAT FLUX INTO TANK, BTU/HR 


VOL * 780. 

MWALL - 975. 

TEMP - -290. 

PRESS « 18. 

PVENT - 30. 

0 - 0 . 

INPUT PARAMETERS FOR FLOW SIMULATION AND COOL-DOWN PIPING* 

CDLENG - LENGTH OF COLL-DOWN PIPING ♦ HEADER 

NNOZ » NUMBER DF NOZZLES 

N07DIA • DIAMETER OF ORIFICE /NOZZLE, FEET 

HEADIA > DIAMETER OF HEADER AND COOL-DOWN PIPING, INCHES 

CONST - 1/ORIFICE CONSTANT. ORIFICE CONSTANT • 0.61 


C 


C 

f 


1 


CDLENG ■ 20. 
NNOZ » 6 
CONST • 1.639 
NOZDIA • 0.5 
HEADIA » 2.0 

GO TO 2 

CALL DATA 02 

MW - 31.9988 
PET ■ 20. 

TET ■ -315. 

OFT - 0. 

VET * 19786. 
TOTAL ■ 6270. 
DIAM - 4. 

NELBOW « 20 
NGATE • 0 
NGLOBE - 2 
NANGLE « 1 
NBUTT ■ 0 
LENGTH - 100. 
METER » 1 
OPIPE « 104400. 
POWER » 4. 

VOL « 300. 

MWALL • 350. 
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TEM° - 60. 

PRESS « 1. 

P VENT - 30. 

0-0. 

CDLENG » 20. 

NNOZ - 6 
CONST ■ 1.639 
NPZDIA ■ 0.25 
HEADIA » 1.0 
C 

2 CONTINUE 

PEAD IN VALUES OF SIMULATION PARAMETERS THAT ARE TO BE DIFFERENT 
FROM THE DEFAULT VALUES. 

READC5.NAME) 

WRITE SOME SIMULATION PARAMETERS TO BE WRITTEN ON DISSPLA PLOTS 
FOR SIMULATION IDENTIFICATION. 

WRITE ( 8»* ) DIAM, POWER, TEMP, PRESS 

WRITE SIMULATION PARAMETERS ON OUTPUT. 

IF(NGAS.E0.2)G0 TO 4449 
WRITE C 6,4400) 

4400 FORMAT (1H1>*HYDR06EN TRANSFER*) 

60 TO 4499 

4449 WRITE(6,4450) 

4450 F0RMAT(1H1,*0XYGEN TRANSFER*) 

4499 WR ITE ( 6, 4500 ) PET»TET»TQTAL,QET 

4500 F0RMAT(1H1,*SUPPLY TANK INITIAL CONDITIONS**/, 5X,*PRE$SURE • * 

♦ ,F10.3,* PSIA*/,5X,*TEMPERATURE • *,F10.3,* F*/,5X, 

♦ *m*$S - *,F10.2,* POUNDS*/, 5X,*HEAT LEAK - *, 

♦ F12.3,* BTU/HR*) 

WRITF (6,4600) PR E SS, TEMP, PVENT, VOL , M VAL L , 0 
4600 FORMATt* RECEIVER TANK INITIAL CONDITIONS **/>5X, 

♦ +PRESSURE - *,F10.3,* PSIA*/,5X»*TEMPERATURE » *,F10.3, 

♦ * F*/,5X,*VENT PRESSURE « *,F10.3,* PSIA*/,5X, 

♦ *TANK VOLUME • *,F10.3,* CUBIC FEET*/»5X, 

*TANK WALL MASS » *,F10.3,* POUNDS*/, 5X,*HEAT LEAK - *, 

♦ F12 . 3, * BTU/HR*) 

WRITE l 6, 4700) DI A M, LENGTH, NELBOW, NG ATE, NGLOBE.N ANGLE, NBUTT, 

♦ METER, OPIPE, POWER 

4700 FORMAT ( * TRANSFER P AR AMETER S » * / , 5 X, *PIPE DIAMETER ■ *, 

♦ FB.3,* INCHES*/, 5X,*LENGTH OF STRAIGHT PIPE - *,F8.3, 

♦ * FEET*/, 5X, 13,* ELBDW(S)*/,5X,I3,* GATE VAL VE ( S )*/ , 5X, 1 3, 

♦ * GLOBE VALVE(S)*/,5X,I3,* ANGLE VAL VE (S )♦/, 5X, I 3, 

♦ * BUTTERFLY VAL VE ( S ) */ ,5 X, 13 ,* FLOW METER ( S )*/, 5X, 

♦ *HE AT LEAK INTO PIPING ■ *,F12.3,* BTU/HR*/,5X, 

♦ *PUMP POWER « *, F5.1,* HP*) 
WRITE(6,4800)C0LENG,NN0Z,N0IDIA,HEADIA 

4800 FORMATt* COOL-DOWN PARAMETERS «*/,5X,*LENGTH OF PIPING ■ * 
1,F8.3,* FEET*/,5X,*NUMBER OF NOZZLES • *, 14, 

1/,5X,*N0ZZLE DIAMETER - *,F8.4,* INCHES*, 

1/5X,*HE ADER DIAMETER • *,F8.4,* INCHES*) 

CALCULATE EQUIVALENT LENGTH OF PIPING. 

LENGTH - LENGTH ♦ FL0AT(NELB0W)*(2.52738*DIAM+. 414286) 

♦ ♦ FL OA T(NGATE)*(. 57083*0 IAM-.0125J+FL0 AT (NGLOBE)* 

♦ (38.0595*DIAM*.10714>*( FLOAT ( NANGL E ) ♦FLOATINBUTT ) /8. ) * 

♦ (14. 3452*DIAM-1. 17857) 

C CONVERT ALL PARAMETERS AND VARIABLES FROM ENGLISH UNITS INTO 


73 



moo o o o o o o o oonno nnr» oo o o o o r> o o o o rt no 


C SI UNITS* PRESSURE - ATM, TEMPERATURE - DEGREES KELVIN* HEAT LEAKS - 
JOULES/SEC# NASS - GRAM-MOLES* PI°E DIAMETER - CM, PIPE 
LENGTH - METERS* PUMP POWER - JCULES/SEC, VOLUME - LITERS* 
TIME - SEC. 

PET - PET/14. 696 
TET • (TET+460. )*5„/9. 

OET - QET*1054. 35/3600. 

VET « VET*28.31? 

TOTAL • T0TAL*454./MW 
DIAM > DIAM*2.54 
LENGTH « LENGTH 4 .3048 
PKETER • FL0AT(METFR)/14.696 
OPIPE » OPIPE*1054.35/3600. 

POWER • POW ER *746 «, 

VOL • VOL*28.317 
MWALL « MWALL*4f 3 »59 
TEMP « (TEMP+460.}*5./9. 

PRESS • PRESS/14.696 
P VENT « PVENT/14.696 
0 » 0*1054.35/3600. 

CDLENG « CDLENG*0 .304 8 
CO I AM « H£A0IA*2. 54 


INITIALIZE THE FLUID PROPERTIES 
OET - FIND D(PET»TET) 

HET « ENTHAL(PET,DET*TET) 


FRACT « TOTAL/VET/DET 
TETSAT « FINOTV(PET) 

ETVAP - F IN DD(PET»TETSAT*1.E-3)*VET*(1. -FRACT) 

OBTAIN OUALITY OF SUPPLY TANK FLUID 

CALL CHECK(HET,PET,TET,ETQ) 

INITIALIZE TIME 

TIME ■ 0. 

INITIALIZE COOL-DOWN NOZZLE PRESSURE DROP 

INITIALIZE FLOWRATE. THE FLOW IS DUE TO PRESSURE DROP BETWEEN THE TANKS 

FC-UESS « FLOtPRESS) 

FL0W1 • PUMPFLO(PRESS*F GUESS) 

INITIALIZE VENT FLOWRATE 

FL0W2 - 0. 

INITIALIZE TOTAL AMOUNT OF MATERIAL VENTED AND AMOUNT OF VAPOR 
PRESENT IN THE RECEIVER TANK. 

X VENT » 0. 

MVAP » 0. 

INITIALIZE COUNTER FOR OUTPUT TO GO TO PLOTTING ROUTINE. 

N PLOT ■ 0 
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C INITIALIZE VENT FLAG TO INDICATE THAT THE RFCFIVER TANK VENT 
C IS CLOSED. 

C 

N VENT - 0 
C 

c initialize programmed stop parameter 
c 

Dl'MP - C. 

c 

c 

C CALL SUBROUTINE THAT MODELS THE COOL-DOWN OF THE RECEIVER TANK UNTIL 
C TWO PHASES ARE PPFSENT IN THE RECFIVFR TANK, 

C 

CALL COOL(MVAP,XVENT*OPRESS) 

NPRT « 0 

c 

C CALCULATE THE ENTHALPY OP THE FLUID ENTERING THE RECEIVER TANK 

C AND CHECK THE OUALITY OF THE FLUID. 

C 

HI ■ HET+OPIPc/FLOWl ♦ P0WER/FLPV1 
CALL CHECK(H1,PPFSS»TEM°»TANK0) 

C 

C SET THE TIME STEP FOR THE SIMULATION. NOTE 1 OTIME • 2C00./FL0W1 
C FOP LH? AND OOO./FLOWl FOR LOT ON NASA-JSC COMPUTER 
C 

DTIMO . VOL /BL0V1 / TOO . 

IFfNGAS .EO. ? ) DT I “E - DTIME/2. 

IFtNGAS ,E0. 1 ) DT I ME « DTIME*2. 

DTTmc . 0TIME+5C. 

C. DTI ME - DTIME/3.1A159 

C INITIALIZE THE AMOUNT OF LIQUID PRESENT IN THE RECEIVER TANK. 

C 

M L 1 0 * 0. 

C 

c 

OLIO » FIND D(PPESS+I.E-3»TEMP) 

DVAP « FIND 0<PPESS-1.E-5»TFMP) 

C 

C LOOP FOR TWO-PHASE SIMULATION. 

I * 0 

1000 CONTINUE 
I ■ I + I 
C 

C OPTAIN VAPOR AND LIOUID ENTHALPY (J/GMQL)ANn INTERNAL ENERGY (J/GMOL) 
C 

H V A P « ENTHAL(PPESS-1.E-B, DVAP, TEMP) 

HlTO • ENTHAKPRFSS+l.E-5, OLIO, TEMP) 

UVAP « HVA°-PPESS*101.3?7/DVAP 
ULIO ■ MLI0-PRESS*101.327/DLI0 
C CALCULATE FUNCTION FOR DELTA P. 

C 

FP « MLIQ*( (l.-DLIQ/nVAP)+DUOPUPRESS)+DLIO+DVDPL(PRESS) 

♦ *(UVAP-l!LIO) )+MVAP*( (l.-DLICZDVAP)*DUDPV(PPESS) 

+ +DLIQ*DVOPV( PR ESS)* (UVAP -ULIO) ) 

F P P « FP + MWALL*CWALL(TEMP)+DTDP<PRESS)* 

♦ (l.-DLIO/DVAP) 

DELTAT » DTDP(PPESS)*DPR£S$ 

C 

C CALCULATE DELTA P 
C 

nPPFS* • ((Hl-UVAP-fDLI0/0VA°*(ULI0-Hl))*FL0Wl + (0)* 

+ (1.-DLI0/DVAP))/FPP 
32 CONTINUE 
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C CALCULATE NEW PRESSURE AND OBTAIN NEW TEMPERATURE. 

C 

PRESS » PRESS «• DPRESS+DTIME 
TEMP « FIND TV (PRESS) 

C 

C RECALCULATE THE FLOWRATE 
FL • FLOWl*MW*. 13216 
FLOW • FL OW1 

FLOW1 - PUMPFLO(PRESStFLOW) 

FL ■ FLOW1+MW*. 13216 
17 CONTINUE 

C CHECK FOR RELIEF VALVE OPENING 
C 

FLOW2 • 0. 

179 CONTINUE 

IF (PRESS +D PR ESS+DTIME .LTo PVENT4.999) GO TO 15 
IF (NVENT .EO. 0 .AND. PR ESS+DPR E SS+DTI ME »LT. PVENT*1.002> 
+ GO TO 15 

C CALCULATE THE VENT FLOWRATE 

FL0W2 • (0*(1.-DIIQ/DVAP)-{UVAP-H1-DLI0/DVAP*(ULI0-H1)) 

♦ *FLOW1)/(DLIO/DVAP*(HVAP~HLIQ) ) 

C 

C CHECK TO SEE IF THE SYSTEM IS VENTING PROPERLY. 

C 

IF (FL0W2 .LT, 0.) GO TO AO 
GO TO 1 A 
AO DPRESS « 0. 

PRESS « PVENT 
C 

C SET FLAG TO INDICATE THAT THE RECEIVER TANK IS VENTING AND 
C CALCULATE TOTAL MOLES VENTED. 

C 

NVENT » 2 

X VENT - XVENT-FL0W2+DTIME 
GO TO 15 

14 FL0W2 - 0. 

NVENT « 0 

C 

15 CONTINUE 
C 

C CALCULATE NEW VAPOR AND LIQUID DENSITIES 
C 

DLIQN - FIND D{PRESSn.E-2»TEMP> 

DVAP • FIND D (PRESS-1. E-5, TEMP) 

IF(DLIQ/DLIQN .GT. lO.JGO TO 231 
DLIC ■ DLIQN 
231 CONTINUE 
C 

C CALCULATE CHANGE IN LIQUID MOLES ANO ACTUAL NUMBER OF MOLES 
C IN LIOUID AND VAPOR PHASES AT T • TIME. 

DMLIQ - -DLI0/DVAP+(FL0Wl*FL0W2)/(l.-DLI0/DVAP)-DLIQ* 

♦ (MLIO*DVDPL(PRF$S)aMVAP*0VDPV( PRESS) )*DPRESS/ 

♦ (l.-DLIO/DVAP) 

MLIO ■ MLIQ+DMLIO*DTIME 

c 

MVAP ■ ( VOL-MLIQ/DlIQ)*DVAP 
C 

C CHECK TO SEE IF THE RECEIVER TANK IS FULL OP OVERFLOWING 
IF(MLIQ4-MVAP + XVENT .LT. 0.95*T0TAL)G0 TO A5 
A6 IF(MLIO+MVAP+XV£NT .l£. TOTAL ) GO TO 48 

47 MLIO - ML I Q-D ML IQ I'D TIME 
MVAP » (VOL -MLIO/DLIQ)* DVAP 

48 DTIME1 « (TOTAL"(XVENT*MLIO+MVAP))/DMLIQ 
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IFtDTIME .17 . DTIHE1 ) 60 TO 45 
49 MLIQ • MLI0*DHLI0*DTIME1 

MVAP • (VOL-MLIO/DLIO)*DVAP 
38 IF (MVAP ,6T. 0.) GO TO 35 

CALL THE SUBROUTINE TO STOP THE PROGRAM. 

NPLOT ■ NPLOT ♦ 1 
X(NPLOT) ■ TIME/60. 

V1CNPL0T) - TENP*9./5.-460. 

Y2(NPL0T> - PRE $S*14. 696 
Y3 (NPLOT) ■ FL0W1+MW*. 13216 
Y4(NPL0T) - M VAP*MW/454. 

Y5 (NPLOT) • MLIOPMW/454. 

Y6(NRL0T) - Y4(NPL0T)+Y5 (NPLOT) 

Y7 ( NPLOT ) • ETO 
YB(NPLOT) » TANKO 
CALL FULL(MLIQ»XVENT) 

CHECK TO SEE IF ALL THE FLUID HAS BEEN TRANSFERRED. 

CALL SUBROUTINE TO STOP PROGRAM. 

35 NPLOT » NPLOT + 1 
X (NPLOT ) ■ TIME/60. 

Yl(NPlOT) • TEMP*9./5,-460. 

Y2 (NPLOT) » PRESS*14.696 
Y3(NPL0T> » FL0W1*MW*. 13216 
Y4 (NPLOT) - MVAPPMW/454. 

Y5(NPL0T) » MLIOPMW/454. 

Y6 (NPLOT) - Y5 ( NPL OT ) ♦ Y4(NPL0T) 

YT(NPLOT) - ETO 

YB(NPLOT) » TANKO 

CALL DONE (MLIO+MVAP, XVENT) 

UPDATE THE SUPPLY TANK PROPERTIES . 

45 ETLIQ ■ TOTAL-MLIO-MVAP-XVENT 
CALL CHECK(HET>PET,TET>ETO) 

CALL SUPPLY (DTIME ) 

CALL CHECK(HET,PET#TET,ETO) 

INCREMENT THE TIME 

TIME • TIME + DTIME 

CHECK TO SEE IF CLOW HAS STOPPED DUE TO RECEIVER TANK BACK PRESSURE 
IN THE ABSENCE OF A PUMP. 

IF (PRESS .GE. PET .AND. POWER .LT. ,05)G0 TO 63 


CHECK THE FLOWRATE TO SEE IF BACK PRESSURE FROM THE RECEIVER HAS 
STOPPED FLOW IN THE PRESECE OF A PUMP. 


IF (ABS(FLOWl) ,LT. 
GO TO 73 
63 CONTINUE 

NPLOT - NPLOT ♦ 1 
X (N PLOT ) • TIME/60. 


,1) GO TO 63 


Yl(NPLOT) 
Y2(NPL0T) 
Y3 (NPLOT) 
Y 4 (NPLOT ) 
Y5(N°L0T) 


TEMP*9./5.-460. 
PRESS*14.696 
F10W14MW*. 13216 
MVAPPMW/ 454. 
MLIQPMW/454. 
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Y6 ( NPLOT) « YMNPL0mY5(NPL0T) 

Y7(HPL0T) - ETO 

Y8CNPL0T) ■ TANKQ 

CALL PFAIL(MVAP*HLIQ»XVENT) 

CALCULATE THE ENTHALPY OF FLUID ENTERING RECEIVER TANK AND OBTAIN 
ITS QUALITY. 

73 CONTINUE 

HI • HET+QPIPE/FL0W1 ♦ P0WER/FLOW1 
CALL CHECKIH1*PRE$S»TEMP, TANKQ) 

IF { DUHP .NE. 1 . )GO TO 585 
NPLOT » NPLOT ♦ 1 
X( NPLOT) • TINE/60. 

Y1 ( NPLOT ) - TEMP**). /5. -460. 

Y2 (NPLOT) « PRES$*14.6<)6 
Y3 ( NPLOT) - FL0V1*MW*.13216 
Y4 { NPLOT) » MVAP*MW/454. 

Y5 { NPLOT) ■ MLIO+MW/454. 

Y6 t NPLOT) » YA(NP|.0T)*Y5(NPL0T) 

Y7 ( NPLOT ) - ETO 

Y8(N°LOT) • TANKQ 

CALL PDATA(NVAP+MLIQ#XVFNT) 

585 CONTINUE 

WRITE RESULTS EVERY CHOSEN TINE INCREMENT. 

NPRT « NPRT ♦ 1 

IF {NPRT .NE. 1) GO TO 20 

NPRT • 0 

NPLOT - NPLOT ♦ 1 
X ( NPLOT ) • TIHE/60. 

Yl(NPLOT) - TEMPT9./5.-460. 

Y2 ( NPLOT) - PRE5S*14.696 
Y3{ NPLOT) ■ FL0W1*MW*. 13216 
YA(NPLOT) « MVAP*MW/454. 

Y5 (N°LOT) • MLI0*MW/454. 

Y6(NPL0T) • YA(NPlOmY5(NPLOT) 

Y7 (NPLOT) ■ ETO 
Y8( NPLOT) - TANKQ 
C 
C 

20 IF(NPLOT ,LE. 600) GO TO 1000 
WR ITE ( 6» 153 ) 

153 FORMAT(/>* ARRAY LENGTH EXCEEDED*) 

CALL PDATA(MVAP+MLIQ»XVENTJ 
END 
C 

C ************ *************************************************************« 

c 

SUBROUTINE COOL(H*XVENT*OPRESS) 

C 

C THIS SUBROUTINE TAKES THE SIMULATION FROM TIME ZERO UNTIL 
C TWO PHASES ARE PRESENT IN THE RECEIVER TANK. 

C 

common/flow/ diam#det#tet# PET>LENGTH,PMETER 
COMMON /PAP AM/ MW,MWALL»Q» VOL »QPIPE»QET» TOTAL »P VENT 
COMMON/SUB/ FL0W1#FL0W2#P«ESS#T,TI«E»HET 
COMMON/NOZZLE/CDIAM,NOZDIA»HEADIA,NNPZ»CONST*CDLENG 
COMMON/PUMP/ POWER# EFF 

COMMON /PLOT/ NPLOT >X»Y1»Y?»Y3#Y4»Y5»Y6»Y7#Y8 

COMMON/ETANK/ETLIO,ETVAP»VET#ETO 

COMMON/PROB/DUMP 
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COMMON/PROP/NGAS 

COMMON/NEG/IJK 

REAL X{ 600) #Y1(600)#Y2(600)»Y3{600)»Y4{600)»Y5(600)#Y6(600) 

♦ » Y7(600 ) # Y8 ( 600) 

REAL MW»MWALL#M#LENGTH#CDLENG#NDZDIA 

INITIALIZE PRINTOUT COUNTER# VENT FLAG AND SIMULATION LOOP COUNTER. 

IJK • 2 
NPP.T - 0 
N VENT ■ 0 
KOUNT ■ 0 
K » 0 

INITIALIZE FLOWRATE BETWEEN TANKS# PUMP EFFICIENCY AND DENSITY OF VAPO 
IN RECEIVER TANK. 

EFF - PEFF(FLOWl) 

DENS • FIND D(PRESS#T) 

CALCULATE THE CHANGE IN VAPOR MOLES. 

ON • FL0W14FL0W2 
DTIME - VOL/FLOW1/8O0. 

IFCNGAS.EQ.2)0TINE-0TIMg/2. 

DTINE - DTIME/7. 

LOOP FOR SIMULATION 

10 KOUNT « KOUNT 4 1 

CALCULATE GRUENSIAN PARAMETER AND THE ENTHALPY OF FLUID ENTERING RECEI 
TANK AND CHECK ITS QUALITY. 

CALL PPHI(DENS#T,PHII) 

GPU - PH 1 1 

HI • HET40PIPE/FL0W1 4 P0WER/FL0W1 
CALL CHECK(H1»PRESS»T#TANK0) 

OBTAIN ENTHALPY OF VAPOR IN RECEIVER TANK (J/GMOL)# HEAT CAPACITY AT 
CONSTANT VOLUME (J/GMOL/DEG K). 

H • ENTHAL(PRESS#DENS»T) 

CVAP - CV(DENS»T) 

CALCULATE CHANGE IN VAPOR DENSITY AND RECEIVER TANK TEMPERATURE 
WITH RESPECT TO TIME. 

DDFNS • DM4DTIME/V0L 

DT-(GRU*T*DM4((H1-H)*FL0W140) /CVAP)/(DENS*V0L4MWALl*CWALl(T) /CVAP) 

4 * DTI ME 

CALCULATE HEAT ABSORBED BY TANK WALL# HEAT OF EXPULSION AND CHANGE 
OF RECEIVER TANK PRESSURE WITH RESPECT TO TIME. 

QWALL - -MWALL*CWALL(T)*DT/DTIME 
THETA « TH(DENS»T) 

0PRESS-GRU*(THETA*DM4(H1-H)*FL0W1404QWALL)*DTIME 
4 /VOL/101.327 

CALCULATE NEW RECEIVER TANK DENSITY# TEMPERATURE AND PRESSURE# AND 
INCREMENT THE SIMULATION TIME. 

DENS ■ DENS 4 ODENS 
T « T 4 OT 
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TIME - TIME ♦ OTIME 
PRESS - PRESS ♦ DPR ESS 
C RECALCULATE THE FLOWRATE 
FLOW « FL0W1 

FLOWI « PUMPFL0(PPESS,FL0W1 
FL • FLOWI* MVP. 13 21 6 
11 CONTINUE 
C 

FLOW? » 0. 

C CHFCK TO SEE IF THE RELIEF VALVE HAS OPENED. 

C 

IF (PRESS .LT. PV£NT*.999! GO TO 25 

IF (NVENT ,E0. 0 .AND. PRESS .LT. PVENT*1.002) GO TO 25 
C 

C CALCULATE THE CHANGE OF MOLES IN RECEIVER TANK WITH VENTING. 

C 

DM « ({H-H11*FLQW1-0-OVALL!/THETA 

CHECK TO SEE IF THE SYSTEM IS VENTING PROPERLY. 

IF ( ( DM-FL0W1 ) .LT. 0.) GO TO 30 
I F ( PRES S .LT. PV£NT*1.002!G0 TO 25 

WR ITE ( 6*2 50 1 

250 FORMAT ( 1H » *THE SYSTEM IS NOT VENTING PROPERLY!** 

♦ * CONTACT A PROGRAMMER*! 

NPLOT • NPLOT ♦ 1 
X(NPLOT) - TIME/60. 

Y1 (NPLOT! « T*9./5.-460. 

Y2 (NPLOT! • PRESS*14.696 
Y3 (NPLOT! • FL0W1*MW*. 13216 
YM NPLOT! • M+MW/454, 

Yf(NPLOT) * 0. 

Y6 (NPLOT! • Y4( NPLOT !*Y5 (NPLOT! 

Y7( NPLOT) - ETO 
Y8 ( N°LOT ) « TANKO 
CALL RD AT A( M, XVFNT > 

30 DPRESS • 0. 

PRESS • PVENT 
FL0W2 ■ OM-FLOWl 

SET FLAG TO SHOW THE RECEIVER TANK IS VENTING AND CALCULATE TOTAL 
AMOUNT VENTED. 

N VENT ■ 2 

XVENT - XVENT-FL0W2*DTIME 
60 TO 26 

25 CONTINUE 
C 

26 DM » FL0W1+FL0W2 
C 

C CALCULATE NEW NUMBER OF MOLES PRESENT IN RECEIVER TANK. 

C 

M « M ♦ DM*DTIME 
C 

C CHECK TO SEE IF ALL FLUID HAS BEEN TRANSFERRED. 

C 

I F ( M+XVENT .LT. 0 . 95*T0T AL 1 GO TO 45 

46 IF(M*XVENT .LE. TOTAL ! GO TO 48 

47 M « M-DM*OT IME 

48 0TIME1 • (TOTAL-(XVENT+M) )/DN 
IF ( DTI ME .LT. DTIME1JG0 TO 4:5 

49 M ■ M+DM*0TIME1 
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CALL DONE (M»XVENT) 

C 

45 CALL CHECK! HE T»PET#TET#ETO) 

UPDATE THE SUPPLY TANK PROPERTIES . 

ETLIQ • TOT AL-H-XVENT 
CALL SUPPLY (DTIME ) 

IF BACK PRESSURE FRDM RECEIVER TANK IN THE ABSENCE OF A PUMP STOPS 
FLOW, CALL SUBROUTINE TO STOP PROGRAM. 

IF (PRESS ,6T. PET .AND. POWER .LT. .05)G0 TO 33 

IF FLOW HAS STOPPED IN THE PRESENCE OF A PUMP DUE TO BACK PRESSURE 
FROM THE RECEIVER# CALL A SUBROUTINE TO STOP PROGRAM. 

IF (ABS(FLOWl) .LT. .1) GO TO 33 
GO TO 23 
33 CONTINUE 

NPLOT • NPLOT ♦ 1 
X(NPLOT) - TIME/60. 

Y1 (NPLOT) • TP9./5.-460. 

Y2 (NPLOT) - PRESS*14.696 
Y3(NPL0T) • FLOWl*MW*. 13216 
Y4(N»L0T) ■ MPMW/454. 

Y5 ( NPL OT) - 0. 

Y6(NPL0T) ■ Y4(NPLOT)+Y5(NPLOT) 

YT(NPLOT) - ETO 
YB(NPLOT) » TANKO 
CALL PF AI L ( M# XVENT ) 

23 CONTINUE 

IF(DUMP .NE. 1 . ) GO TO 485 
NPLOT ■ NPLOT ♦ 1 
X (NPLOT ) ■ TIME/60. 

Y1 ( NPLOT) • T*9. /5.-460. 

Y2 ( NPLOT ) - PRESS*14.696 
Y3( NPLOT) » FLOWl+MW*. 13216 
Y4 (NPLOT) « MPMW/454. 

Y5CNPLOT) • 0. 

Y6 (NPLOT ) ■ Y4(NPL0T)+Y5(NPL0T) 

Y? ( N p LOT ) - ETO 
YB(NPLOT) • TANKO 
CALL PDATAfM, XVENT) 

485 CONTINUE 

NPRT ■ NPRT ♦ 1 

WRITE SIMULATION VALUES INTO PLOT ARRAYS EVERY CHOSEN ITERATION, 

IF (NPRT .NE. 5) GO TO 35 
NPRT • 0 
K - 0 

NPLOT ■ NPLOT ♦ 1 
X (NPLOT) - TIME/60. 

Y1 (NPLOT) • T*9./5.-460. 

Y2(NPL0T) • P PE SS*14. 69 6 
Y3 (NPLOT ) ■ FL0W1*MW*. 13216 
Y4 ( NPL OT) • MPMW/454. 

Y5 ( NPLOT) • 0. 

Y 6 (NPLOT) - Y4 ( NPLOT ) 

YT(NPLOT) - ETQ 
YB(NPLOT) « TANKO 

IF RECEIVER TANK TEMPERATURE IS GREATER THAN CRITICAL TEMPERATURE# 
DON'T BOTHER CHECKING FOR TWO PHASES IN RECEIVER TANK. 
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IF(NGA$.EQ.2) GO TO 19 
35 IF (T . GE . 40. > GO TO 20 
GO TO 36 

19 IF(T.GE.157.)G0 TO 20 
C 

C OBTAIN SATURATION TEMPERATURE AT RECEIVER TANK PRESSURE# USE THIS TO 
C GET SATURATION DENSITY} COMPARE SIMULATION AND SATURATION DENSITIES 

C TO DETERMINE WHETHER OR NOT ANY LIQUID IS PRESENT IN THE RECEIVER TANI 

C 

36 TSAT ■ FIND TV(PRESS) 

TTPRES - PRESS-1. E-5 

DS AT - FIND D(TTPRIS»TSAT) 

DCALC • M/VOL 
IF (T .LT. TSAT) T • TSAT 
IF (DCALC .IE. DSAT ) GO TO 20 
C 

C IF TWO PHASES ARE PRESENT, BACK UP ONE ITERATION AND RETURN TO THE 
C main PROGRAM. SHUT OFF COOL-DOWN NOZZLES# RECALCULATE FLOWRATE. 

C 

M « M - DM*DTIME 
DENS » DENS-DOENS 
TIME ■ TIME-DTIME 
PRESS - PRESS-DPRESS 
T « FIND TV ( PRESS ) 

C 

CONST « 0. 

CDLENG • 0. 

IJK « 1 

FGUESS « FL0W1 * 5. 

IFCNGAS .EQ. 2)FGUESS • FGUESS*4. 

FL0W1 ■ PUNPFLO(PRESS#F GUESS) 

C 

TIM « TIME/60. 

WRITE(6#58)TIM#K0UNT 

58 FORMAT ( /# * COOLDOWN TIME - *#l : 12.4#* KCUNT ■ *#I6) 

RETURN 

20 I F (NPLOT .L E. 600)60 TO 10 
WR ITE (6,151 ) 

151 FORMAT ( * ARRAY LENGTH EXCEEDED*) 

CALL PDATA( M# XVENT) 

END 

C 

C ********************************** *************************************** 

c 

subroutine PPHKDD# TT#PHII) 

c 

C SUBROUTINE TO CALCULATE THE GRUENS IAN PARAMETER (PHII) AS A 
C FUNCTION OF DENSITY (6M0L/L) AND TEMPERATURE (DEG K). 

C 

D « OD 
T • TT 

CALL DPDT (DT» D» T > 

CAP - C V(D»T) 

PHII - DT/D/CAPU01.325 
RETURN 
END 
C 

C******************** **************************************** ************* 

C 

FUNCTION TH ( D ENS# TEMP ) 

C 

C CALCULATE THE HEAT OF EXPULSION (J/GMOL) AS A FUNCTION OF 
C DENSITY (GMOL/L) AND TEMPERATURE (DEG K). 
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C 


CALL DPDT(DT»DENS*TEMP) 

CALL DPDD(DD*DENS*TEMP) 

TH • CP(DENS,TEMP)*DENS*OD/DT 

RETURN 

END 

****** ** ******************************** ************ ********************** 
FUNCTION OVDPV(PRESS) 

COMMON/ PROS/ DUMP 

CALCULATE THE PARTIAL DERIVATIVE OF SPECIFIC VOLUME WITH RESPECT TO 
PRESSURE FOR A SATURATED VAPOR. 

PVAP ■ PRESS-1. E-5 
TEMP . FINDTV( PRESS) 

V • l./FIND D ( PVAP* TEMP ) 

SET AN INITIAL DELTA P 

DP » .05 

IF< (PRESS-DP) • LT. 1.E-AJ0P - DP/2. 

MOUNT • 0 

25 TEMP - FIND TV(PRESS-DP) 

VLOV » l./FIND D(PVAP-DP»TEHP) 

TEMP - FIND TV( PRESS+DP ) 

VHIGH « l./FIND 0(PVAP*DP*TEMP) 

DVDPV « (VHIGH-VLOW) /DP/2, 

SEE IF THE DERIVATIVE WILL CHANGE BY MAKING DELTA P SMALLER 

IF C ABS ( ( VHIGH-2 •* V*VLQW) /( VHIGH-V) ) .LT. l.E-3) RETURN 
KOUNT ■ MOUNT ♦ 1 

IF THE SOLUTION HAS NOT CONVERGED, MAKE DELTA P SMALLER. 

DP ■ OP /2. 

IF THE DERIVATIVE HAS NOT CONVERGED IN 100 ITERATIONS, WRITE AN 
ERROR MESSAGE (THIS FUNCTION USUALLY CONVERGES WITHIN 2-3 ITERATIONS) 

IF (MOUNT ,LE. 100) GO TO 25 
WRITE (6*125) 

125 FORMATdH * *KQUNT EXCEEDS 100 FOR DVDPV. STOP PROGRAM*) 

DUMP - 1. 

RETURN 

END 


******* 4 **** **************************************** ********************* 


FUNCTION DVDPL(PRESS) 

COMMON/PROB/DUMP 

CALCULATE THE PARTIAL DERIVATIVE OF SPECIFIC VOLUME WITH RESPECT TO 
PRESSURE FOR SATURATED LIQUID, INPUT IS PRESSURE (ATM), OUTPUT 
DVDPL (L/GMOL/ATM) 

PLIO « PR ES S+l.E-5 
TEMP . FIND TV(PRESS) 

V ■ l./FIND D(PLIO*TEMP) 

DP - .05 
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IF( (PRESS-DP) .LT. X.E-*)DP « DP/2. 

KOUNT - 0 

25 TEMP • FIND TV(PRESS-DP) 

VLOW • 1. /FIND D(PLIO-DP,TEMP> 

TEMP ■ FIND TV(PRESS*DP) 

VHIGH ■ l./FIND D(PLIQ*DP,TEMP) 

DVDPL - { VHIGH-VLOW) /DP/2. 

SEE IF THE PARTIAL DERIVATIVE HAS CONVERGED. 

IF (AB$((VHIGH-2.*V+VL0V>/(VHIGH-V)> .LE. l.E-3) RETURN 
KOUNT ■ KOUNT* 1 

IF THE SOLUTION HAS NOT CONVERGED. MAKE DELTA P SMALLER. 

DP « DP / 2. 

THE DERIVATIVE SHOULD CONVERGE IN 100 ITERATIONS (USUALLY, 
2-3 ITERATIONS ARE SUFFICIENT?. 

IF (KOUNT .LE. 100) GO TO 25 
VR ITE (6,125) 

125 FORMAT < 1H ,*KOUNT EXCEEDS 100 FOR DVDPL. STOP PROGRAM*) 
DUMP • x. 

RETURN 

END 


********** ************************************************ ************** 

FUNCTION DUDPV(PRESS) 

C0MM0N/PR08/DUMP 

CALCULATE THE PARTIAL DERIVATIVE OF INTERNAL ENERC-Y WITH RESPECT 
TO PRESSURE FOR A SATURATED VAPOR. INPUT IS PRESSURE (ATM). 

TEMP • FINDTV ( PPESS ) 

DENS • FIND D(PRESS-1.E-5,TEMP> 

SINCE INTERNAL ENERGY IS NOT RETURNED BY THE THERMODYNAMIC PROPERTIES 
PACKAGE, U ■ H - PV 

U • ENTHAL(PRESS-l.E-5, DENS, TI:HP)-PRESS*101. 327/DENS 
DP - .05 

I F ( (PRESS-DP) .LT. l.E-4)DP ■ DP/2. 

KOUNT « 0 

25 TEMP « FIND TV (PRESS -DP) 

DENS ■ FIND D ( PRESS-0 P-1. E-5, TEMP) 

ULOW - ENTHALt PRESS-DP-1. E-5, DENS, TEMP )-(PRESS-DP)*101. 327/DENS 
TEMP - FINDTV(PRE$S+DP) 

DENS ■ FIND D(PRESS*DP-1. E-5, TEMP) 

UHIGH « ENTHAL ( PR F SS+DP-l. E-5, DENS, TEMP )-{ PRESS+DP) *101. 327/DENS 
DUOPV • ( UHIGH-ULPW) /DP/2 • 

SEE IF THE DERIVATIVE HAS CONVERGED. 

IF ( AB S ( (UHI GH-2 ,*U*ULOW) /(UHIGH-U)).LT, l.E-3) RETURN 
KOUNT ■ KOUNT + 1 

IF THE SOLUTION HAS NOT CONVERGED, MAKE DELTA P SMALLER. 

DP ■ OP /2. 

IF (KOUNT .LE. 100) GO TO 25 
WR ITE (6,125) 
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125 FORMAT ( IN »*K0UNT EXCEEDS 100 FOR DUDPV. STOP PROGRAM*) 

DUMP ■ 1. 

RETURN 

END 

C 

c** **** ****** ********* * * ***************** ********************************** 

c 

FUNCTION DUDPL (PRESS) 

C 

COMMON/ PROB/DUMP 

CALCULATE THE PARTIAL DERIVATIVE OF INTERNAL ENERGY WITH RESPECT TO 
PRESSURE FOR SATURATED LIOUIO. INPUT IS PRESSURE (ATM), PARTIAL 
DERIVATVE OUTPUT ( J /GMOL / ATM) . 

TEMP « FI NDTV < PRESS ) 

DENS • FIND D(PRESS+1.E-5»TEMP) 

U • ENTHAL(PRESS+l.E-5»DENS»TEMP)-PRESS*101. 327/DENS 
DP • .05 

IF( (PRESS-OP) .LT.1.E-A)DP ■ DP/2. 

KOUNT ■ 0 

25 TEMP ■ FI NDTV (PRESS- DP) 

DENS - FIND 0 ( PPESS-DP*1 .E-5.TEMP) 

UlOW * ENTHAL { PRE SS-DP + 1 • E-5»DENS» TEMP )-(PRESS-DP)*101, 327 /DENS 
TEMP - FINDTV(PPESS*DP) 

DENS - FIND D(PRESS*DP*1.E-5»TEMP) 

UHIGH • ENTHAL ( PR ESS +DP+1.E-5. DENS, TEMP )-( PRESS + DP )*101.327/ DENS 
DUDPL ■ (UHIGH-UL0W)/DP/2. 

CHECK ACCURACY OF PARTIAL DERIVATIVE. 

IF(ABS( (UHI GH— 2 .+U+-ULOW ) / { UHI GH-U ) ) »LT • l.E-3) RETURN 
KOUNT • KOUNT ♦ 1 

IF THE SOLUTION HAS NOT CONVERGED. MAKE DELTA P SMALLER. 

DP « DP /2. 

CHECK TO SEE IF THE SOLUTION IS CONVERGING (CONVERGENCE USUALLY 
OCCURS IN 2-3 ITERATIONS). 

IF (KOUNT .LE. 100) GO TO 25 
WRITE (6.125 ) 

125 FORMAT ( 1H .*KOUNT EXCEEDS 100 FOR DUDPL. STOP PR06R AM* ) 

DUMP • 1. 

RETURN 
END 


******«****************4***t***********************#*******************t*> 

FUNCTION DTDP ( PRESS ) 

C CALCULATE THE PARTIAL DERIVATIVE OF TEMPERTURE WITH 
C RESPECT TO PRESS FOR A SATURATED LIOUID 
C 

COMMON/ PR08/0UMP 
C 

TEMP « FINDTV(PRESS) 

DP - 0.05 

IF(PRESS-DP .LT. l.E-4)DP«DP/2. 

KOUNT » 0 

25 THIGH ■ F INOTV( PR E SS+DP ) 

TLOW » FINDTV(PRESS-DP) 

DTDP • (THIGH-TLOW) /DP/2, 

C 
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C CHECK CONVERGENCE 
C 

I F C ABS C (THIGH-2. +TEMP+TLOW) /( THIGH-TEMP ) ) ,L£. l.E-3) 

♦ RETURN 
KOUNT ■ KOUNT+1 
C 

C REDUCE DP 
C 

DP - DP/2. 

C 

C 

IF (KOUNT ,LE. 100) GO TO 25 
WRITE (6*125 ) 

125 F0RMAT(1H **KOUNT EXCEEDS 100 FOR DTDP. STOP PROGRAM*) 

DUMP « 1. 

RETURN 

END 

C **************************************************** ********************** 
C 

FUNCTION CVAll (TEMP ) 

C 

C CALCULATE THE HEAT CAPACITY OF THE TANK WALL (J/G/DEG K) 

C AS A FUNCTION OF TEMPERATURE ( OEG K). WALL MATERIAL IS ALUMINUM 
C FOR LH2 AND INCONEL FOR LOX. 

C 

COMMON/PROP/NGAS 
IF(NSAS»EQ»2) GO TO 1 

CW ALL • EXP(7. 29296-121. 912 /TEMP+347. 91 8/ (TEMP**2 ) > /1000. 

GO TO 2 

1 CWALL • EXP(6. 482-106. 3/TEMP+321.9/(TEMP**2))/1000. 

2 RETURN 

END 

C 

C**** ******** ****** ********************** ************** ******************** 

c 

FUNCTION FLO( PRESS ) 

C 

C CALCULATE THE PUMPLESS FLOWRATE (GMOL/SEC) AS A FUNCTION OF 
C SUPPLY AND RECEIVER TANK PRESSURES (ATM). FOR THIS FUNCTION* THE 
C SUPPLY TANK PRESSURE MUST BE GREATER THAN THE RECEIVER TANK PRESSURE. 

C 

COMMON /FLOW/ D X AM, DET>TET, PET* LENGTH* PME TER 
CDMM0N/N0ZZLE/CDIAM*N0ZDXA*HEADIA*NNC7*C0NST*CDLENG 
COMMON/PARAM/ MW, MWALL » 0, VOL* QPI PE* QET* TOTAL, P VENT 
C0KM0N/PR08/DUMP 
REAL LENGTH, MH»MWALL»NOZDIA 
PI - 3.14159 
C 

C INITIALIZE THE FRICTION FACTOR. 

C 

FLINE • 0.005 
FCOOL - 0.005 
PNETER • (PET-PRESS)*0.05 
KOUNT • 0 
C 

C CALCULATE FLOWRATE* OBTAIN A NEW FANNING FRICTION FACTOR* THEN 

C RE-CALCULATE THE FLOWRATE. THIS CONTINUES UNTIL THE FRICTION 

C FACTOR CHANGES S IGN IF IC ANTLV 

C 

35 CONTINUE 
C 

C CALCULATE CONSTANTS FOR INDIVIDUAL PRESSURE DROPS 
C 

C MAIN TRANSFER LINE 



c 

CLINE - (3.15ei5*MW*LENGTH)/<DET*PI**2*DIAM**5) 

C 

C COOLDOWN LINE 
C 

CCOOL ■ (3.15815*MW*CDLENG)/{DET*PI*+2*CDIAM**3) 

C 

C COOLDOWN NOZZLES 
C 

A - CONST * 6.7162E-2 / (NOZDIA*+2) 

B • 1.941E-3 ♦ MW / OET 
C - 1. - ((N0ZDIA/HEADIA)**4) 

CNOZ - ( { A**2 ) * B * C > / FLOAT(NNOZ) * 4.725E-4 
C 

C CALCULATE FLOWRATE 
C 

FLO - $ QRT( ( PET-PR ESS-PMETER ) / C CLINE*FL INE+CCOQL*FCOOL*CNOZ ) ) 

C 

C CALCULATE NEW FRICTION FACTORS 
C 

FLNEW • FRICT<FLO*TET»DET,DIAM) 

FCNEW • FRICT(FLO*TET»OET,CDIAM) 

C 

C CONVERGENCE CHECK 
C 

IFf ABSIFLNEW-FLINE ) /FLINE . LT. O.C1 . AND. A BS ( FCNE W-FCOOL ) /FCOOL 
♦ .LT.O.ODGO TO 80 
GO TO 61 
60 CONTINUE 

C. WRITE (6#77)FLINE#FC00L 

77 FORMAT (* FLINE » +»E12.5»* FCOOL « **E12.5> 

RETURN 

81 CONTINUE 
C ITERATE 
C 

FLINE - FLNEW 

FCOOL - FCNEW 

KOUNT • KOUNT+1 

IF (KOUNT .LT. 100)G0 TO 35 

WPITE(6>100 ) 

100 FpRfAT ( 1H t AKOUNT EXCEEDS 100 FOR FLO. STOP PROGRAM* ) 

101 DUMP - 1. 

RETURN 

END 

C 

£************************************************************************** 

c 

FUNCTION FRICT(FLO*T*D*DIAM) 

C 

C CALCULATE THE FRICTION FACTOR AS A FUNCTION OF PIPE ROUGHNESS AND 
C REYNOLDS NUMBER. 

C 

CONNON/PARAH/ NW*MWALL»0»V0L»0PIPE*0ET,T0TAL»PVENT 
CONMON/PROP/NGAS 
RFAL NW.MWALL 
PI - 3.14159 
C 

C CALCULATE THE DIMENSIONLESS REYNOLDS NUMBER 

C 

IFtNGAS.E0.2) GO TO 1 

REY • 4.*FL0PMW/(1.8367E-4*DIAM*PI) 

GO TO 2 

1 REY » 4.*FL0*MW/(VISC<D#T)*1.E-6*DIAH*PI) 

C 
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C THE PIPE ROUGHNESS IS ASSUMED TO BE 0.00015 INCHES 
C 

2 EPSILON • .00015*30.48 

FRICT - (-4.*AL0G10(EPSIL0N/0IAM/3.7-5 .C2/REY* AL0G101 EPSI LON 
♦ /DIAM/3.7 * 14.5/REY)>)**{-2 > 

RETURN 

END 

C 

C ********************************** ********************************** *******' 

c 

SUBROUTINE CHECK(H1»PRESS, TEMP, OUAL) 

CHECK THE OUALITY OF A STREAK OR TANK GIVEN ENTHALPY (J/GMOL), 

PRESSURE (ATM) AND TEMPERATURE (DEG K). 

TSAT ■ FIND TV ( PRESS ) 

DS AT - FIND D(PRESS*1.E-5,TSAT) 

HCALC • ENTHAL(PRE$S+1.E-5,D$AT,TSAT> 

IF (HCALC .LT. HI) GO TO 10 

IF ACTUAL ENTHALPY IS LESS THAN LIOUID SATURATION ENTHALPY, THE 
LIQUID IS SUBCOOLED AND THERE IS NO VAPOR PRESENT (QUAL - 0). 

OUAL - 0. 

RFTURN 

10 DVAP - FIND D(PRESS-l.E-5»TSATil 

HVAP • ENTHALtPRESS-l.E-5, DVAP, TSAT) 

CALCULATE THE QUALITY USING ACTUAL, LIOUID AND VAPOR ENTHALPIES. 

OUAL »1 00 •*(HCALC-Hl)/{ HCALC- HVAP) 

RETURN 
END 


* ********* ************************************** ****** ********* 


FUNCTION PUMPFL0(PRESS»PFL01) 

CALCULATE THE FLOWRATE BETWEEN TANKS (GMOL/SEC) IN THE PRESENCE OF 
A PUMP (J/S). THIS ROUTINE REQUIRES AN INITIAL GUESS FOR THE 
FLOWRATE. 

CO«MON/FLOW/ DIAM,DET,TET, PET, LENGTH, PMETER 
COMMON/NOZZLE/CDIAN,NOZDIA,HEADIA,NNOZ,CONST,CDLENG 
COMMON /PAR AM/ MW, M WALL , 0, VOL, QPI PE , GET, TOTAL , P VENT 
COMMON/PUMP/ POWER, EFF 
COMMON/PROB/DUMP 
COMMON/NEC-/IJK 
REAL MW, MWALL, LENGTH, NOZDIA 
KOUNT ■ 0 
DP • PET-PRESS 
PMETER « DP*0.05 
PI • 3.14159 
25 EFF « PEFF(PFLOl) 

CALCULATE A THEORETICAL PRESSURE DROP (OR RISE) BETWEEN TANKS. 

NOZZLE PRESSURE DROP 

A • CONST * 6.7162E -2 / (N0ZDIA**2) 

8 • 1.941E-3 * MW / DET 
C • 1. - <(N0ZDIA/HEADIA)**4) 

DPNOZ » ( ( P FL01**2) * (A+*2) *8*0/ FLOAT (NNOZ) * 4.725E-4 
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C COOLDOWN LINE 
C 

FCOOL • FRICT(PFLOl.TET.DET.CDIAM) 

D PC POL* (PFL 01**2 )*(3»15815*MW*C OLE NG*FC COL )/(DET*PI**2*CDIAM**5) 
MAIN TRANSFER LINE 

FLINE « FRICT(PFLG1»TET»DET»DIAM) 

DPLINE • (PFL01**2)*(3.15815*MW*LENGTH*FLINE) / (DET*P1**2*DIAM**5) 
PUMP 

OPPUMP » (0.009868*POVEP*EFF*DET>/(PFL01) 

CALCULATE PRESSURE DROP 

DPCALC1 - ♦ PHETER ♦ DPNOZ ♦ DPC00L ♦ DPLINE - OPPUMP 

INCREMENT FOR INTER ATION 

MOUNT - MOUNT ♦ 1 
IF (MOUNT .GT. 1) GO TO 10 
PFLO « PFL01 
PFL01 - PFL0in.5 
DPCALC • DPCALC1 
GO TO 25 

CALCULATE A NEW FLOWRATE USING A OUASI-NEWTON METHOD. 

NOMENCLATURE I 

PUMPFLO - FLOWRATE (GMOL/S) AT N*1 ITERATION 
PFLOl - FLDWARATE AT N ITERATION 
PFLO - FLOWATE AT N-l ITERATION 

DPCALC1 - CALCULATED PRESSURE DROP (ATM) AT N ITERATION 
DPCALC - CALCULATED PRESSURE DROP AT N-l ITERATION 
DP - ACTUAL PRESSURE DRO° (OR RISE) BETWEEN TANKS. 

10 PUMPFLO - PFL01-( DP-DPC ALC1)* (PFL01-PFL0)/(DPCALC-DPCALC1) 

CHECK TO SEE IF THE THEORETICAL AND ACTUAL TANK PRESSURE 
DIFFERENCES ARE WITHIN 0.01? OF EACH OTHER. 


CONVERGENCE CHECK 

I F ( ASS C (0PCALC1-DP)/DP) .LT. l.E-A) GO TO 180 
GO TO 170 
180 CONTINUE 
160 CONTINUE 
RETURN 

170 CONTINUE 

ITERATE THE FLOWRATE. CALCULATE PRESSURE DROP. 

DPCALC • DPCAIC1 
PFLO » PFLOl 
PFLOl ■ PUMPFLO 
C 

IF (MOUNT .LT. 100) GO TO 25 
WP ITE ( 6.100 ) 

100 FORNATUH ,*KOUNT EXCEEDS 100 FOR PUMPFLO. STOP PROGRAM*) 
DUMP . 1. 

RETURN 

END 
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C 

C** ********** ********** ************************ ************** **************** 

c 

FUNCTION PEFF ( FL0W1 ! 

CALCULATE PUMP EFFICIENCY OF A CENTAUR PUMP 

COMMON/FLOW/ DI AM, DET, TET, PET, L ENGTH, PMFTER 
COMMON/PROP/NGAS 
REAL LENGTH 

CHANGE UNITS ON FLOWRATE FROM GNOL/S TO GALLONS/MIN. 

GPM • FL0W1*15.8502/DET 
CALCULATE PUMP EFFICIENCY 
IF ( NGAS . E 0« 2 ) GO TO 1 

PEFF-(.39213*.09662*6PM*2.7797E-5*(GPM**2)- 

♦ 8. 1723 E- 8* (GPM*+3)+2.5008E-ll7(6PM**6))/100. 

GO TO 2 

1 PEFF » (.32031+. 48261*GPM-1.6223E-3*(GPM**2)* 

♦ 3. 9836 E-6* ( GPM**3 ) -4.885 3E-9* ( GPM**4 ) )/100. 

2 CONTINUE 

RETURN 
END 

*************************************************************************** 


SUBROUTINE PF AI L ( M , X VENT ) 

THIS SUBROUTINE FINISHES OFF THE PROGRAM. PFAIL IS USED WHEN 
FLOW BETWEEN TANKS STOPS DUE TO BACK PRESSURE FROM THE RECEIVER 
TANK. 

COMMON /SUB/ FLOWl,FLOW2, PR ESS, TEMP, TIME, HET 
COMMON /PAR AM/ MW, MWALL , 0, VOL , OPI PE , QET, TOTAL , P VENT 
COMMON /FLOW/ D I AM , DET, TET, PE T, LENGTH , PMETE R 
COMMON/P LOT/ NPL0T,X,Y1»Y2»Y3»Y4,Y5»Y6»Y7»Y8 
COMMON/PROP/NGAS 
DIMENSION OPLDT(IB) 

REAL X(600) ,Y1(600),Y2(600) »Y3(600), Y4(600),Y5 (600)»Y6(600) 

♦ ,Y7(600),Y8(600) 

REAL M, MW, MWALL, LENGTH 

CONVERT TIME FROM SECONDS TO MINUTES, PRESSURE FROM ATM TO PSI, 
TEMPERATURE FROM DEG K TO DEG F AND QUANTITY FROM GMOL TO LB. 

XTIME ■ TIME/60. 

X PRESS • PRESS*14.696 
PET • PET*14« 696 
XTEMP » TEMP79./5.-460. 

TET • TET79./5.-460. 

X M » M*MW/454. 

XOUT » XVENT*MW/454. 

WRITE A MESSAGE STATING WHY THE SIMULATION STOPPED AND WHAT THE 
CONDITIONS WERE AT THAT TIME. 

WRITE (6, 100) XTIME, XTEMP, XPRESS,XM, XOUT 
100 FORMATt//* AT *,F7.3,* MINUTES FROM START OF PROPELLANT* 

+ ,* TRANSFER*/,* FLOW HAS STOPPED DUE TO BACK PRESSURE FROM* 

♦ ,* THE RECEIVER TANK.*, //»*AT THIS TIME THE RECEIVER TANK*, 

♦ * IS AT ♦, F7 .2,* DEG F, *,F7»2,* PSIA*/, 
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♦ * AND CONTAINS *,F10.2,* POUNDS *, 

♦ //,* DURING THE SIMULATION *,F8.2,* POUNDS HERE VENTED*) 

C 

VRITE!6,101)TET,PET 

101 FORMAT!/* SUPPLY TANK IS AT *F9.A,* DEG F, *F7.2,* PSIA*) 

SKIP OVER OTHER MESSAGES AND GO TO INITIALIZATION OF PLOT DATA 
AND PARAMETERS. 

GO TO 10 
ENTRY FULL 

THIS ENTRY IS USED WHEN THE RECEIVER TANK IS FULL OF LIOUID. 

CONVERT TIME, TEMPERATURE, PRESSURE AND OUANTITY IN RECEIVER TANK TO 
ENGLISH UNITS. 

IFCM . L T. TOTAL .AND. XVENT .EO. 0.)M ■ TOTAL 
I F ( M+XVENT .NE. TOTAL ) M - TOT AL-XVENT 
XTIME » TIME/60. 

XTEM» • TEMP*9./5.-A60. 

TET • TET+9./5.-A60. 

XPRESS ■ PRE$S*1A.696 
PET » PET*1A.696 
XM « M*MH /A 56 • 

XOUT - XVENT*MW/A5A • 

WRITE OUT A MESSAGE STATING WHY THE SIMULATION STOPPED AND HHAT THE 
CONDITIONS HERE AT THAT TINE. 

WRITE <6,110) XTIME, XTEMP, XPRESS , XM, XOUT 

110 FORMAT!//,* AT *,F7.3,* MINUTES FROM START OF PROPELLANT*, 

♦ * TRANSFER,*/,* THE RECEIVER TANK IS FULL.*//,* AT THIS TIME * 

♦ *THE RECEIVER TANK IS AT *»F7.2»* DEG F, *»F7.2»* PSIA,*/, 

♦ * AND CONTAINS *,F10.2,*P0UNDS*, 

♦ //,* DURING THE SIMULATION *,F8.2,* POUNDS HERE VENTEO*) 
WRIT£(6»111 )T£T,P£T 

111 FORMAT!/* SUPPLY TANK IS AT *F9.A,* DEG F, *F7.2»* PSIA*) 

SKIP OVER UNWANTED MESSAGE. 

GO TO 10 
ENTRY DONE 

THIS ENTRY IS USED WHEN THE SUPPLY TANK IS EMPTY. 

I F ( M ,LT, TOTAL .AND. XVENT .EQ. 0.)N ■ TOTAL 
IF ( M+XVENT .NE. TOTAL ) M « TOTAL-XVENT 
CONVERT TO ENGLISH UNITS. 

XTIME • TIME/60. 

XTEMP ■ TEMP*9./5.-A60. 

TET - TET*9./5.-A60. 

XPRESS • PR ESS*1A • 696 
PET - P ET*1 A. 696 
XM ■ M*MW/A5A» 

XOUT • XVENT*MW/A5 A. 

WRITE AN APPROPRIATE MESSAGE AND STOPPING CONDITIONS. 

WRITE (6, 120) XTIME, XTEMP, XPRESS, XM, XOUT 
120 FORMAT!/,* AT *,F7.3,* MINUTES FROM START OF PROPELLANT*, 

♦ ♦ TRANSFER*/,* THE SUPPLY TANK IS EMPTY*,/,* AT THIS TIME THE* 

♦ * RECEIVER TANK IS AT *,F7.2,* OFG F, *,F7.2,* PSIA,*/, 

♦ * AND CONTAINS *, F10 .2, *POUNDS*» 
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♦ /,* DURING THE SIMULATION *»F8.2, 

♦ * POUNDS WERE VENTEO*) 

C 

WRITE(6,121)TET,PET 

121 FORHATt /* SUPPLY TANK IS AT *F9.4»* DEG F, *F7.2,* PS I A* ) 

C IDENTIFY THE PROPELLANT FOR THE PLOTTING ROUTINE. 

C 

GO TO 1C 
ENTRY PDATA 
C 

C THIS ENTRY IS USED WHEN A PROGRAMMED STOP IS EXECUTED 
C 

C CONVERT TO ENGLISH UNITS 
C 

XTIME ■ TIME/60. 

XPRESS • PR ES$*14 .696 
PET ■ PET*14.696 
XTEMP - TEMP*9./5.-460. 

TET • TET*9./5.-460. 

XM > M*MW/454. 

XOUT • XVENT*MW/454. 

C 

C WRITF A MESSAGE STATING THE CONDITIONS WHEN PR06RAM IS STOPPED 
C 

WRITE(6»140)XTIME,XTEMP,XPR£SS,XM,X0UT 

140 FORMAT ( / / »* AT 4, FT. 3,* MINUTES FROM START OF PROPELLANT*, 

♦ * TRANSFER*/,* A STOP IS REACHEO*, //* AT THIS TIME THE* 

♦ * RECEIVER TANK IS AT *,F7.2»* DEG F, *,F7.2»* PSIA,*/, 

♦ * AND CONTAINS *, F10. 2,* POUNDS *, 

♦ //,* DURING THE SIMULATION *,F8.2, 

♦ * POUNDS WERE VENTED*) 

C 

WRITE( 6,141 )TET, PET 

141 FORMAT!//* SUPPLY TANK IS AT *F9.4,* DEG F, *F7.2»* PSIA*) 

10 CONTINUE 

I F (NG AS .E 0. 2 ) GO TO 11 
C 

WP ITE ( 8,190 ) 

190 FORMAT(*HYDROGEN*) 

GO TO 12 

11 WR ITE (8,191 ) 

191 FORMAT(*OXYGEN*) 

12 CONTINUE 
C 

C WRITE THE NUMBER OF PLOT POINTS FOR USE BY THE PLOT ROUTINE. 
C 

WP ITE ( 8, 200 ) NPLOT 
200 FORMAT { 13 ) 

C 

C DPL0T(4) ARE BOUNDS FOR Yl, 0PL0T<5) AND OPLOT<6) ARE BOUNDS 
C SET LOWER AND UPPER BOUNDARIES FOR PLOTS 
C 

OPLOT(l) • 0. 

0PL0TC2 ) ■ 0. 

DO 708 1*1,4 

I F ( TIME /60. .LE. 10,)0 p LQT(2) - 10, 

I F ( TI HE /60» .LE. 5.)0PL0T(2) • 5. 

IF (TIME/60. .GT. OPLOT (2 ) )0P LOT ( 2) • 0PL0T<2> ♦ 50. 

708 CONTINUE 

OPLOT (3 ) » -500. 

IFCNGAS.EQ.2)0PL0T(3> ■ -400. 

OPLOT (4 ) * 0. 

IF(NGAS.E0.2)0PL0T(4) - 100. 

DO 55 J » 1,3 
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IF(Y1(1> .GT. OPLOT (4 ) ) 0PL0T(4) • 0PL0T(4 >+100. 

55 CONTINUE 

OPLOT ( 5 ) - 0. 

0PL0TC6) • 50. 

DO 57 I-1»NPL0T 

IF({Y2(I).GT.0PL0T(6)).AND.(Y2(I).LT.(0PL0T{6)+50.)))0PL0T(6) • 

+ OPLOT { 6) 4 50. 

57 CONTINUE 
56 0PL0TC7) ■ 0. 

OPLOT < 8 ) - 2000. 

IF(NGAS.EQ.2)0PL0T(8) • 10000. 

OPLOT (9 ) • 0. 

OPLOmOJ - 250. 

TOTAL - TOTAL/454. *MW 
DO 1000 1-5000,65000,5000 

IF( TOTAL, LE.FL0AT(I))0PL0T( 10) • FL0AT<I>*0.10 
IF ( TOTAL. LE.FL0AT(I))G0 TO 1001 

1000 CONTINUE 

1001 CONTINUE 
OPLOT ( 11 ) - 0. 

OPLOTf 12) ■ OPLOT ( 10>*10. 

0PL0T113) • 0. 

OPLOT (14) • OPLOT { 12 ) 

OPL OT( 1 5 ) - 0. 

OPLOT (16) • 100. 

OPLOT (17) - 0. 

0PL0T(18) - 100. 

WRITE PLOT PARAMETERS FOR PLOTTING ROUTINE. 

WRITE (6, ♦) (OPLOT(II), 11-1,18) 

WRITE(8,*)(0PL0T{II),II-1,18) 

WRITF ARRAYS CONTAINING SIMULATION VARIABLES. 

WRITE(6,249)(X(I),Y1(I),I-1,2) 

249 FORMAT ( 2F15 ,5 ) 

WRITE (8,250) ( X( I ) , Y1 ( I ), Y2 { I ), Y3 ( I ), Y4 ( I ) , Y5 ( I ), Y6 ( I ) , 

4 Y4CI), Y8(I),I»1,NPL0T) 

250 F0RMAT(F8.3,T10,F8.3,T20,F8.5,T30,F8.2,T40,F8.1,T50,F8.1,T60,F8.1 
4 ,T70,F8.3,T80,F8.3) 

STOP 

END 

****** ******* ******************************************** ******* **#*****■ 


SUBROUTINE SUPPLY ( DTIME ) 

COMMON /SUB/ FL0W1,FL0W2,PRESS,TEMP,TIME,HET 
COMMON /FLOW/ DIAM,DET,TET, PET, LENGTH, PMETER 
COMMON /PAR AM/ MW, MWALL , 0, VOL, QP I PE , OET, TOTAL, P VENT 
C0MM0N/PR08/DUMP 

COMMON /E TANK /ETLIO,ETVAP, VET, ETO 
REAL MW, MWALL, LENGTH 

THIS SUBROUTINE MONITORS THE PROPERTIES OF THE SUPPLY TANK. 
INITIALIZE THE SUBROUTINE 

ENTHALPY AND HEAT CAPACITY OF LIOUIO IN SUPPLY TANK. 

HL « HET 

CSUBP • CP(DET,TET> 

C FRACTION OF TANK VOLUME OCCUPIED BY LIQUID AND AMOUNT OF 
C HEAT TO EACH PHASE. 
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FR ACT - ETLIQ/VET /OFT 
DOL « OET*FRACT 
OQV • OET-DOL 
C 

C SPECIFIC VOLUME OF LIQUID . 

C 

VFTLIQ ■ l./DET 
C 

C SATURATION TEMP AT GIVEN PRESSURE* AND THE VAPOR DENSITY. 

C 

TETSAT - FINDTV(PET) 

DFTVAP « FINDD<PET,TETSAT+l.E-3> 

C 

C SPECIFIC VOLUME OF VAPOR. 

C 

VETVAP ■ l./DETVAP 
C 
C 

C GRUNEISEN PARAMETER FOR vapor AND LIQUID. 

C 

CALL PPHI(DETVAP*TETSAT«-1.E-3»PHII) 

GRUVAP « PHII 

CALL PPHI(DET»TET*PHII) 

GRULIQ - PHII 
C 

C HEAT OF EXPULSION FOR VAPOR AND LIQUID 
C 

THETAV - TH(DETVAP,TET$AT+l.E-3) 

THETAL « TH(DET»TIET) 

C 

C CALCULATE THE CHANGE IN THE SUPPLY TANK PRESS WITH TIME 
C 

DPET ■ ({-FLOW1+VETVAP/VFTL IQ/THE TAV*DQV*1./THETAL 
1 *DOL) /ETLIO/(VETLIQ/THETAL/GRULIO+ETVAP* 

? VETVAP**2/ETLIQ/VETHQ/THETAV/GRUVAP)/101.327)*DTIM£ 

C 

C CALCULATE THE CHANGE IN THE SUPPLY TANK LIQUID TEMPERATURE 
C WITH TIME 
C 

DTET « ({OOL/ETLIQ+VETLIQ*ALPHAL(PET»T£T* DET)*DPET*101 .327 
1 /DTIME )/CSUPP)*DTIME 
C 

C CALCULATE THE CHANGE IN THE SUPPLY TANK LIQUID ENTHALPY WITH 
C TIME 
C 

DHL - (DQL/ETLI 0+VETLIQ*D PET /DTIME *101.327) *DTIME 
C 

C INCREMENT THE PRESSURE* LIQUID TEMPERATURE AND ENTHALPY 
C 

PET » PET+DPET 
TET - TET+DTET 
HL • HL+DHL 

C OBTAIN NEW LIQUID DENSITY AND ENTHALPY FROM FLUIDSPACK 

C OBTAIN THE SUPPLY TANK LIQUID QUALITY 

C 

501 DET - FINDD(PET*TET) 

HET « ENTHAHPET*DET»TET> 

CALL CHECK(HET,PET,TET,ET0) 

C 

C IF OUALITY IS GREATER THAN 7ER0* CAVITATION COULD OCCUR. STOP 
C TRANSFER IMMEDIATELY. 

IFt ETQ.LT .1 .E-6)G0 TO 99 
WR ITF ( 6*225 ) ETO 
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225 F0RMAT(1H0»*ET LIO BOILING. STOP TR AN$FER**5X* *ETO ■ *E12.5) 

DUMP « 1. 

RETURN 

COMPARE MODEL AND FLUIDSPACK ENTHALPIES. IF THE VALUES ARE WIDELY 
DISPARATE* THE ET MODEL IS IN ERROR. STOP THE PROGRAM. 

OTHERWISE USE THE MODEL VALUE. 

99 IF(A8S( (HET-HL)/HL) .GT. O.10O)GO TO 100 
HET « HL 
RETURN 

IOC WRITE (6*255 )HET*HL 

255 FORMAT { 1H0* * ENTHALPY DISPARITY. HET « *E12.5** HL ■ *E12.5) 

DUMP - 1. 

RETURN 

END 

**** ** ******************************** ********* ***************** ****** 
FUNCTION ALPHAL(PET*TET*DET) 

COMN0N/PROB/OUMP 

THIS FUNCTION CALCULATE THE BULK EXPANSIVITY FOR THE LIQUID 


SET THE TEMPERATURE STEP AND THE LOOP COUNTER 

DT « 0.05 
KOUNT ■ 0 

FIND VALUES SLIGHTLY ABOVE AND BELOW THE ACTUAL TEMPERATURE 
FOR THE DENSITY 

25 DHIGH - FINDD(PET,TET-DT) 

DLOW - FINOD(PET,TET+DT> 

CALCULATE THE BULK EXPANSIVITY 

ALPHAL • -TET/DET*(DHIGH-0L0W)/DT/2, 

I F ( APS { (DHIGH-2«*DET+DL0W)/ ( DHIGH—DET ) ) .LT. l.E-3)RETURN 
KOUNT « KOUNT ♦ 1 

IF THE VALUES ARE NOT CONVERGING* MAKE THE TEMPERATURE STEP SMALLER 
DT • DT/2. 

IF(K0UNT.LE.100)G0 TO 25 
WR ITE (6*125 ) 

125 FORMAT (1H1»*K0UNT EXCEEDS 100 FOR ALPHAL. STOP PROGRAM*) 

DUMP - 1. 

RETURN 

END 
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.PROC ,DEMO . 

.* 

.* DEMONSTRATION PROGRAM USING DISSPLA 9.0. 

EVICT, LGO, PLOT. 

FTN , I =PRO G , L=FAI L , R=3 , PMD . 

IFE,FILE(GRAFING, . NOT. AS) , JO. 

GET , GRAF ING/UN=G RAF . 

ENDIF , JO . 

GET, TAPE2=PSDATA. 

BEGIN, DISPLA9 ,GRAFING . 

RETURN, PROG, LGO. 

REPLACE, PLOT. 

RENAME , Z=PLOT . 

SKIPR,Z . 

COPYBR , Z , OUTPUT , 1 . 

REVERT. 

EXIT. 

REVERT, ABORT. 

REPLACE, FAIL. 

.DATA, PROG 

PROGRAM DEMO (INPUT, TAPE 2 ,OUTPUT=/1000 ,PL0T=/1000) 

REAL LMIN, LMAX 

DIMENSION X(400) ,Y1(400) ,Y2(400) ,Y3(400) ,Y4(400) ,Y5(400 ) , 

+ Y6(400) ,Y7 (400) ,Y8(400) 

DIMENSION XAXIS(4) , HEAD ( 3 ) ,IPKRAY(300) ,L1(3) ,L2(3) ,L3(3) ,L4(3) 
DIMENSION NAME(2) ,IPR2(100) 

DATA XAXIS / 1 OHTIME ( MINU , 10HTES ) $ ,2*10H 

READ ( 2 , * ) DIAM, POWER, TEMP, PRESS 
READ(2 ,80) NATE ( 1 ) , NAME ( 2 ) 

80 FORMAT (2A4) 

READ ( 2 , * ) N 

READ ( 2 , * ) TMl , TM2 , TMIN , TMAX , PMIN , PMAX , FMIN , FMAX , VMIN , 

+ VMAX , LMIN , LMAX , DUM1 , DUM2 , DUM3 , DUM4 , DUM5 , DUM6 
DTIME = (TM2-TM1 ) /10 . 

DTEMP = ( TMAX -TMIN) /5 . 

DPRESS = (PMAX-PMIN) /5 . 

DFLOW = (FMAX- FMIN) /5 . 

DVAP = (VMAX-VMIN)/5. 

DLIQ = (LMAX -LMIN) /5 . 

DO 10 1=1, N 

10 READ(2 ,*) X(I),Y1(I),Y2(I),Y3(I),Y4(I),Y5(I),Y6(I),Y7(I),Y8(I) 
CALL ID ( " LOUIE, $”,100) 

C CAL HP7221 

CALL FR80(0 .0 ) 

CAL TK4010(120) 

CALL HWSPEC (4HPLOT , 4HFILE) 

CALL NOBRDR 
CALL PHYSOR( 1 .25,6 . ) 

CALL AREA2D(6. ,4.) 

ENCODE (30, 100, HE A) NAME(l) , NAME (2) 
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100 FORMAT(2A4," PROPELLANT TRANSFER$ ") 

CALL HEADIN(HEAD, 100, 1.25,1) 

CALL XNAME(XAXIS.IOO) 

CALL YNAME( "TEMPERATURE (DEG F)$",100) 
CALL YAXANG(0 . ) 

ENC0DE(30 , 15 ,L1 ) DIAM 
15 FORMAT ("DIAMETER: ”,F4.1," IN.$") 

ENCODE (30, 20 ,L2) POWER 
20 FORMAT ("PUMP: ",F4.1," HP$") 
ENCODE(30,25,L3) TEMP 
25 FORMAT ( "TANK TEMP: ",F6.1," DEG F$”) 
ENC0DE(30 ,30 ,L4) PRESS 
30 F0RMAT( "TANK PRESS: ",F4.1," PSIA$”) 

CALL LINES (Ll ,IPKRAY,1) 

CALL LINES(L2 , IPKRAY, 2) 

CALL LINES (L3, IPKRAY, 3) 

CALL LINES(L4, IPKRAY, 4) 

CALL LSTORY ( IPKRAY , 4 , 3 . 5 , 3 . 0 ) 

CALL GRAF( TM1 ,DTIME , TM2 , TMIN , DTEMP , TMAX) 
CALL CURVE(X,Yl,N,0) 

CALL ENDGR(O) 

CALL PHYS0R(1 .25 , 1 . ) 

CALL AREA2D(6. ,4.) 

CALL XNAME(XAXIS , 100) 

CALL YNAME( "PRESSURE (PSIA)$",100) 

CALL YAXANG(0 . ) 

CALL GRAF ( TM1 , DT IHE , TM2 , PMIN , DPRESS , PMAX) 
CALL CURVE(X,Y2,N,0) 

CALL ENDGR(O) 

CALL ENDPL(O) 

CALL NOBRDR 
CALL YAXANG(0 . ) 

CALL PHYS0R(1 .25,6.) 

CALL AREA2D(6. ,4.) 

CALL HEADIN(HEAD ,100,1.25,2) 

CALL HEADIN( "(PAGE 2)$" ,100 ,1. ,2) 

CALL XNAME(XAXIS , 100) 

CALL YNAME( "POUNDS IN VAPOR PHASE? ",100) 
CALL GRAF ( TM1 , DT IME , TM2 , VMIN , DVAP , VMAX ) 
CALL CURVE(X,Y4 ,N,0) 

CALL ENDGR(O) 

CALL YAXANG(0 . ) 

CALL PHYSOR( 1.25,1.) 

CALL AREA2D ( 6 . , 4 . ) 

CALL XNAME (XAXIS.lOO) 

CALL YNAME ("POUNDS IN LIQUID PHASE? ”,100) 
CALL GRAF ( TM1 , DTIME , TM2 , LMIN , DL IQ , LMAX ) 
CALL CURVE(X,Y5,N,0) 

CALL ENDGR(O) 

CALL ENDPL(O) 

CALL NOBRDR 
CALL YAXANG(0 . ) 

CALL PHYS0R(1 .25,6.) 

CALL AREA2D ( 6 . , 4 . ) 
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CALL HEADIN( HEAD, 100,1.25, 2) 

CALL HEAD IN( "(PAGE 3)$ " ,100 ,1 . ,2) 

CALL XNAME(XAXIS ,100) 

CALL YNAME( "TOTAL POUNDS TRANSFERRED$ ” , 100) 
CALL GRAF(TM1,DTIME,TM2,LMIN,DLIQ,LMAX) 

CALL CURVE(X,Y6 ,N,0) 

CALL ENDGR(O) 

CALL YAXANG(O-) 

CALL PHYSOR( 1.25,1.) 

CALL AREA2D(6. ,4.) 

CALL XNAME (XAXIS.100) 

CALL YNAME ("FLOWRATE (POUNDS/MIN)$ " , 100) 
CALL GRAF(TM1 , DTIME , T M2 , FMI N , D FLO W , FMAX ) 

CALL CURVE (X,Y3,N,0) 

CALL ENDPL(O) 

CALL PHYS0R(1 .5,6.) 

CALL AREA2D(6. ,4.) 

CALL HEADIN(HEAD , 100 ,1.25,2) 

CALL HEADIN( ” (PAGE 4) $" ,100 , 1 . ,2 ) 

CALL XNAME(XAXIS , 100 ) 

CALL YNAME (1H ,1) 

CALL YNONUM 

CALL LINES( "SUPPLY$" ,IPK2 , 1) 

CALL LINES( "RECEIVER INLET$ " , IPK2 , 2 ) 

CALL GRAF(TM1, DTIME, TM2, -10. ,120. ,110.) 

CALL RLVEC(TMl- .2 ,0 . ,TMl+.2 ,0 . ,0) 

CALL RLVEC(TMl-.2 , 100. ,TM1+. 2 , 100 . ,0) 

XI = XINVRS (XMESS( "100% LIQ. $”,100) ,0.) 

X2 = XINVRS (XMESS( "100% VAP. $",100),0.) 

CALL RLMESS( "100% LIQ. $" ,100 ,TML- .2-X1 ,0 . ) 
CALL RLMESS( "100% VAP. $" ,100,TMl-.2-X2,100.) 
CALL CURVE (X,Y7,N, 10) 

CALL CURVE(X,Y8 ,N, 10) 

CALL LEGEND ( IPK2 ,2,3. , 3 . ) 

CALL ENDPL(O) 

CALL DONEPL 

STOP 

END 
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PROGRAM PLOT ( INPUT, TAPE6 ,OUTPUT=TAPE6 ,TAPE8=INPUT) 

REAL TIME(400) ,TEMP(400) ,PRESS(400) ,FLOWl(400) ,MVAP(400) ,MLIQ(400) 
REAL LLB, LUB 

READ(8,*) DUM1 , DUM2 , DUM3 , DUM4 
READ(8 ,90) NAME 
90 FORMAT(AIO) 

READ(8 ,100) N 
100 FORMAT (13) 

READ ( 8 , * ) TLB , TUB , TEMPLB , TEMPUB , PLB , PUB , FLB , FUB , VLB , VUB , LLB , LUB , 

+ TOTLB , TOTUB , Q1LB , Q1 UB , Q2LB , Q2UB 
DO 10 I ■ 1,N 

10 READ(8 ,110) TIME(I) ,TEMP(I) ,PRESS(I) ,FL0W1(I) ,MVAP(I) ,MLIQ(I) 

110 F0RMAT(F8 .3 ,T10 ,F8 .3 ,T20 ,F8 .5 ,T30 ,F8 .2 ,T40 ,F8 . 1 ,T50 ,F8 . 1) 

ILOG = 0 
ISIZE = 0 
NOUT - 0 
WRITE (6 ,120) 

120 FORMAT( 1H1 , *ORBITER TANK TEMPERATURE (DEGREES F) VS.* 

+ ,* TIME (MINUTES)*) 

CALL PPC ( TEMP , T IME , N , ILOG , I S I ZE , NOUT , TEMPLB , TEMPUB , TLB , TUB ) 

WRITE(6 , 130) 

130 FORMAT (1H1,*TANK PRESSURE (PSIA) VS. TIME (MINUTES)*) 

CALL PPC ( PRES S , TIME , N , ILOG , I S IZE , NOUT , PLB , PUB , TLB , TUB ) 

WRITE(6 ,140) 

140 F0RMAT(1H1 ,* FLOWRATE (LB/MIN) VS. TIME (MINUTES)*) 

CALL P PC ( FL0W1 , TIME, N , ILOG , ISIZE, NOUT , FLB , FUB , TLB , TUB ) 

WRITE(6 , 150) 

150 F0RMAT(1H1,*P0UNDS IN VAPOR PHASE VS. TIME (MINUTES)*) 

CALL PPC (MVAP , TIME , N , ILOG , ISIZE ,N0UT , VLB , VUB , TLB , TUB) 

WRITE(6 , 160) 

160 FORMAT (1H1,*P0UNDS IN LIQUID PHASE VS. TIME (MINUTES)*) 

CALL PPC(MLIQ, TIME, N, ILOG, ISIZE, NOUT, LLB, LUB, TLB, TUB) 

STOP 

END 
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PSDATA 


4 . 4 . 60 . 1 . 

OXYGEN 
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5 . - 400 . 100 . 

o 

• 

Ln 

O 

• 

O 

• 

10000 . 0 . 

1000 . 0 . 10000 . 0 . 10000 . 0 . 100 . 

0 . 100 

.001 

48.699 

2.75189 

3398.11 

3.1 

0.0 

3.1 

0.000 

.002 

37.590 

4.42648 

3398.11 

6.2 

0.0 

6.2 

0.000 

.003 

26.672 

6.02568 

3398.11 

9.3 

0.0 

9.3 

0.000 

.004 

15.938 

7.55134 

3398.11 

12.4 

0.0 

12.4 

0.000 

.005 

5.384 

9.00520 

3398.11 

15.5 

0.0 

15.5 

0.000 

.005 

- 4.998 

10.38895 

3398.11 

18.6 

0.0 

18.6 

0.000 

.006 

- 15.211 

11.70422 

3398.11 

21.7 

0.0 

21.7 

0.000 

.007 

- 25.260 

12.95254 

3398.11 

24.8 

0.0 

24.8 

0.000 

.008 

- 35.151 

14.13541 

3398.11 

27.9 

0.0 

27.9 

0.000 

.009 

- 44.889 

15.25425 

3398.11 

31.0 

0.0 

31.0 

0.000 

.010 

- 54.477 

16.31042 

3398.11 

34.1 

0,0 

34.1 

0.000 

.011 

- 63.921 

17.30522 

3398.11 

37.2 

0.0 

37.2 

0.000 

.012 

- 73.224 

18.23991 

3398.11 

40.3 

0.0 

40.3 

0.000 

.013 

- 82.391 

19.11569 

3398.11 

43.4 

0.0 

43.4 

0.000 

.014 

- 91.425 

19.93369 

3398.11 

46.5 

0.0 

46.5 

0.000 

.015 

- 100.331 

20.69502 

3398.11 

49.6 

0.0 

49.6 

0.000 

.016 

- 109.113 

21.40072 

3398.11 

52.8 

0.0 

52.8 

0.000 

.016 

- 117.773 

22.05177 

3398.11 

55.9 

0.0 

55.9 

0.000 

.017 

- 126.316 

22.64913 

3398.11 

59.0 

0.0 

59.0 

0.000 

.018 

- 134.745 

23.19370 

3398.11 

62.1 

0.0 

62.1 

0.000 

.019 

- 143.064 

23.68633 

3398.11 

65.2 

0.0 

65.2 

0.000 

.020 

- 151.276 

24.12782 

3398.11 

68.3 

0.0 

68.3 

0.000 

.021 

- 159.383 

24.51894 

3398.11 

71.4 

0.0 

71.4 

0.000 

.022 

- 167.390 

24.86040 

3398.11 

74.5 

0.0 

74.5 

0.000 

.023 

- 175.299 

25.15287 

3398.11 

77.6 

0.0 

77.6 

0.000 

.024 

- 183.114 

25.39698 

3398.11 

80.7 

0.0 

80.7 

0.000 

.025 

- 190.836 

25.59330 

3398.11 

83.8 

0.0 

83.8 

0.000 

.026 

- 198.469 

25.74237 

3398.11 

86.9 

0.0 

86.9 

0.000 

.026 

- 206.015 

25.84467 

3398.11 

90.0 

0.0 

90.0 

0.000 

.027 

- 213.476 

25.90064 

3398.11 

93.1 

0.0 

93.1 

0.000 

.028 

- 220.856 

25.91069 

3398.11 

96.2 

0.0 

96.2 

0.000 

.029 

- 228.155 

25.87516 

3398.11 

99.3 

0.0 

99.3 

0.000 

.030 

- 235.376 

25.79435 

3398.11 

102.4 

0.0 

102.4 

0.000 

.031 

- 242.520 

25.66852 

3398.11 

105.5 

0.0 

105.5 

0.000 

.032 

- 249.589 

25.49788 

3398.11 

108.6 

0.0 

108.6 

0.000 

.033 

- 256.583 

25.28260 

3398.11 

111.7 

0.0 

111.7 

0.000 

.034 

- 263.502 

25.02281 

3398.11 

114.8 

0.0 

114.8 

0.000 

.035 

- 270.347 

24.71859 

3398.11 

117.9 

0.0 

117.9 

0.000 

.036 

- 277.115 

24.37000 

3398.11 

121.0 

0.0 

121.0 

0.000 

.037 

- 283.805 

23.97707 

3398.11 

124.1 

0.0 

124.1 

0.000 

.037 

- 290.413 

23.53982 

3398.11 

127.2 

0.0 

127.2 

0.000 

.038 

- 295.604 

23.06581 

3398.11 

130.3 

0.0 

130.3 

0.000 

.145 

- 291.025 

21.25714 

2628.10 

116.0 

293.3 

409.3 

0.000 

.253 

- 292.181 

19.97737 

2628.10 

108.0 

583.2 

691.3 

0.000 

.360 

- 293.033 

19.07144 

2628.10 

102.1 

871.1 

973.2 

0.000 

.467 

- 293.684 

18.40162 

2628.10 

97.4 

1157.8 

1255.2 

0.000 

.575 

- 294.195 

17.88860 

2628.10 

93.5 

1443.6 

1537.1 

0.000 


8.601 

5.954 

4.107 

2.671 

1.493 

.494 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 
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.682 - 294.606 17.48423 

.789 - 294.942 17.15792 
.896 - 295.223 16.88937 

1.004 - 295.461 16.66470 

1.111 - 295.664 16.47408 
1.218 - 295.840 16.31039 
1.325 - 295.994 16.16835 
1.433 - 296.130 16.04397 
1.540 - 296.251 15.93416 

1.647 - 296.358 15.83652 
1.755 - 296.455 15.74914 
1.862 - 296.543 15.67050 
1.969 - 296.622 15.59936 
2.076 - 296.694 15.53469 
2.184 - 296.761 15.47565 

2.291 - 296.822 15.42154 

2.355 - 296.867 15.38141 


2628.10 

90.3 

1728.8 

2628.10 

87.4 

2013.5 

2628.10 

84.9 

2298.0 

2628.10 

82.6 

2582.2 

2628.10 

80.5 

2866.3 

2628.10 

78.5 

3150.1 

2628.10 

76.7 

3433.9 

2628.10 

74.9 

3717.6 

2628.10 

73.3 

4001.2 

2628.10 

71.7 

4284.7 

2628.10 

70.1 

4568.2 

2628.10 

68.6 

4851.6 

2628.10 

67.2 

5135.0 

2628.10 

65.7 

5418.4 

2628.10 

64.3 

5701.7 

2628.10 

63.0 

5985.0 

2628.10 

61.9 

6211.6 


1819.0 

0.000 

0.000 

2101.0 

0.000 

0.000 

2382.9 

0.000 

0.000 

2664.8 

0.000 

0.000 

2946.8 

0.000 

0.000 

3228.7 

0.000 

0.000 

3510.6 

0.000 

0.000 

3792.5 

0.000 

0.000 

4074.5 

0.000 

0.000 

4356.4 

0.000 

0.000 

4638.3 

0.000 

0.000 

4920.3 

0.000 

0.000 

5202.2 

0.000 

0.000 

5484.1 

0.000 

0.000 

5766.0 

0.000 

0.000 

6048.0 

0.000 

0.000 

6273.5 

0.000 

0.000 
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NOPLOT 


/JOB 
/NOSEQ 
COMBO (T500) 

/READ, ACCOUNT 
HEADING. OLOUIE 

GET, COMBO, BOXB,TAPE5=DATAIN,TAPE7=FUEL. 
FTN,I=C0MB0,0PT=2,L=TRACE,R=3,PMD. 

LDSET ( PRESET =ZERO ) 

LOAD.BOXB. 

LGO. 

SKIP, PASS. 

EXIT. 

REWIND, TRACE. 

COPYEI .TRACE, TAPE6. 

ENDIF.PASS. 

REWIND, TAPE8. 

COPYEI, TAPE8.PSDATA. 

REPLACE, PSDATA. 

DAYFILE.TAPE6. 

REWIND, TAPE6. 

COPYEI ,TAPE6,PS0UT. 

REPLACE, TAPE6.PSOUT. 

DAYFILE.DAYF. 

REPLACE, DAYF. 

/EOI 


102 



LIBPLOT 


/JOB 

/NOSEQ 

FPL0T(T500) 

/READ .ACCOUNT 
GET.PSDATA.FPLOT . 

ATTACH , STAR4 /UN=CAMLIB , NA . 
LIBRARY, STAR4. 

FTN , I=FPLOT , B=B ,L=0 ,R=3 . 
LOAD.B. 

EXE CUTE,,,, PSDATA. 

/READ, SPY 


SUBCOM 


/JOB 

/NOSEQ 

COMBOCT500) 

/RE AD, ACCOUNT 
HEADING. OLOUIE 

GET , COMBO , BOXB , TAPE 5 =D AT AI N , TAPE7=FUEL . 
FTN ,I=COMBO, OPT*2 , L=TRACE , R=3 , PMD . 
LDSET(PRESET=ZERO) 

LOAD, BOXB. 

LGO. 

SKIP, PASS. 

EXIT. 

REWIND, TRACE. 

COPYEI .TRACE ,TAPE6 . 

ENDIF.PASS. 

REWIND, TAPE8. 

COPYE I , TAPE8 , P SDATA . 

REPLACE, PSD AT A. 

REWIND, TAPE8. 

ATTACH , STAR4/UN=CAMLIB ,NA. 

LIBRARY , STAR4 . 

GET,FPLOT. 

FIN , I=PLOT , B=B ,L=LIST,R=3 . 

LOAD , B . 

EXE CUTE,, TAPE8. 

SKIP, OK. 

EXIT. 

REWIND, LIST. 

COPYEI, LIST, TAPE6. 

ENDIF.OK. 

DAYFILE ,TAPE6 . 

REWIND, TAPE6 . 

COPYEI ,TAPE6 ,PSOUT . 

REPLACE , TAPE6 , P SOUT . 

DAYFILE, DAYF. 

REPLACE ,DAYF . 

/EOI 
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DATAIN 


$NAME 

QPIPE = 1000000., 

$ 

FUEL 


$FUEL 

NGAS=2, 

$ 
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Oxygen Default Parameters 


Supply tank initial conditions: 

Pressure = 

Temperature = 

Oxygen Mass = 

Heat Leak = 

Receiver tank initial conditions: 

Pressure = 

Temperature = 

Vent Pressure = 

Tank Volume = 

Tank Wall Mass = 

Heat Leak = 

Transfer Parameters: 

Pipe Diameter = 

Length of Straight Pipe = 

20 Elbow(s) 

0 Gate Valve(s) 

2 Globe Valve(s) 

1 Angle Valve(s) 

0 Butterfly Valve(s) 

1 Flow Meter (s) 

Heat Leak Into Piping = 

Pump Power = 

Cool-Down Time = 

Cool-Down Parameters: 

Length of Piping = 

Header Diameter = 

Nozzle Diameter = 

Number of Nozzles = 


20.000 psia 
-315.000 F 

6270.00 Pounds 
0.000 Btu/hr 


1.000 psia 

60.000 F 

30.000 psia 

300.000 Cubic Feet 

350.000 Pounds 
0.000 Btu/hr 


4.000 Inches 

100.000 Feet 


104400.000 Btu/hr 

4.0 HP 

1.000 Minutes 


20.000 Feet 
0.5000 Inches 
0.1250 Inches 
6 


At 1.92 minutes from start of propellant transfer the supply tank is empty. 

At this time the, receiver tank is at -310.17 Deg F, 6.65 psia, and contains 
6270.00 pounds of oxygen. 

During the simulation 0.00 pounds were vented. 
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OXYGEN 


PROPELLENT IRRNSFER * 
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FLOWRRTE (POUNDS/MIN) *10 TOTRL POUNDS TRRNSFERRED 


OXYGEN PROPELLENT TRRNSEER 

( PAGE 3) 
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OXYGEN PROPELLENT TRRNSFER 

(PAGE 4) 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 


TIME (MINUTES) 
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PRESSURE (PSIRJ TEMPERATURE (DEG F) 


HYDROGEN PROPELLENT TRANSFER 
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TOTAL POUNDS TRANSFERRED 


HYDROGEN PROPELLENT TRRNSFER 

(PAGE 3) 



HYDROGEN PROPELLRNT TRRNSPER 

(PAGE 4) 
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NOMENCLATURE 


The units used in the computer model and derivation of equations are given 
first and are followed by SI units, which are enclosed in parentheses. 


bp 

C 

C 

C' 

C, 


P II 


w 

D 

dm 

f 

F(P) 

h 

L 

L e 

M 

m 

w 


m 

N 

P 

AP 

AP 

AP 


Re 


cooldown 


meter 

Q 

^total 

T 


brake horsepower for a pump, HP, (W) 

conversion factor used in Equations (28) and (29) 

constants used in Equation (17) 

constant used in Equation (30) 

constants used in Equation (39) 

constant pressure heat capacity, Btu/(lb*°F), (J/(gmol*K)) 

constant volume heat capacity, Btu/(lb-°F), (J/(gmol*K)) 

constant pressure heat capacity of receiver tank wall, 
Btu/(lb;°F), (J/(gmol-K)) 

pipe diameter, inches, (cm) 

differential mass or flowrate of fluid, Ib/min, (gmol/s) 

Fanning friction function factor, dimensionless 

term defined in Equation (lib) 

specific enthalpy, Btu/lb, (vl/gmol) 

pipe length, ft, (m) 

equivalent length of a pipefitting, ft, (m) 
molecular weight, Ib/mol, (g/gmol) 
tank wall mass, lb, (g) 
propellant quantity, lb, (moles) 

Reynolds number, dimensionless 
pressure, atm, (MPa) 

pressure drop between supply and receiver tanks, atm, (MPa) 

pressure drop in the cooldown piping, atm, (MPa) 

pressure drop across flowmeter, atm, (MPa) 

heat flux in the receiver tank, Btu/hr, (W) 

total heat flux in the system, Btu/hr, (W) 

temperature, °F, (K) 
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u 

u 

u b 

V 

v 

w 

x 

X 

a 

e 

n 

e 

M 

P 

<l> 

Subscripts 

2 

v 

1 

2 

cool 

i 

mai n 
o 

RT 

s 

sat 

w 


total internal energy, Btu, (J) 
specific internal energy, Btu/lb, (J/gmol) 
bulk velocity, ft/s, (m/s) 
tank volume, cu.ft, (£) 
specific volume, cu. ft/lb, (£/gmol) 
work, Btu/hr, (J/s) 
quality of fluid, dimensionless 
dummy variable used in Equation (15) 
bulk thermal expansivity, dimensionless 

pipe roughness, ft, (m) 
pump efficiency, dimensionless 
heat of expulsion, -p^ 9 ^ 


apj,, Btu/lb, (J/gmol) 
viscosity, lb/(ft*s), (gmol/(m*s)) 
density, lb/ft 3 , (gmo1/£) 

Gruneisen parameter, i(|E) p = f(g) s _ d1mens1 


onless 


1 iquid 
vapor 

main fill line into receiver tank 

receiver tank vent line 

cool down piping 

either stream 1 or 2 

main transfer piping 

orifice in cooldown piping 

receiver tank 

supply tank 

saturated conditions 

tank wall 
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